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PREFACE 


The following compilation of documents includes a list of the 66 attendees, a copy of the 
viewgraphs presented, and a summary of the discussions held after each session at the 1996 
Coolant Flow Management Workshop held at the Ohio Aerospace Institute, adjacent to the 
NASA Lewis Research Center, Cleveland, Ohio on December 12-13, 1996. The workshop was 
organized by H. Joseph Gladden and Steven A. Hippensteele of NASA Lewis Research Center. 
Participants in this workshop included Coolant Flow Management team members from the 
NASA Lewis, their support service contractors, the turbine engine companies, and the 
universities. 

The NASA Lewis Coolant Flow Management Team began their process in response to 
the 1993 Turbomachinery Peer Review committee suggestion and after visiting several aircraft 
engine manufacturers and discussing what was of interest to industry (ref. NASA TM 107161, 
“Research Strategy for Modeling the Complexities of Turbine Heat Transfer” by Robert J. 
Simoneau, March, 1996). Three areas of turbine study developed: (1) details of the internal 
passages, (2) computational film cooling capabilities, and (3) the effects of heat transfer on both 
sides. Several research projects, contracts, and grants have resulted. Also mentioned in this 
reference, this workshop was planned to keep the team focused and to assess the progress of the 
work. 


The purpose of the workshop was to assemble the team members, along with others who 
work in gas turbine cooling research, to discuss needed research and recommend approaches that 
can be incorporated into the Center’s Coolant Flow Management Program. The workshop was 
divided into three sessions: (I) Internal Coolant Passage Presentations, (II) Film Cooling 
Presentations, and (III) Coolant Flow Integration & Optimization. Following each session there 
was a group discussion period. 


Steven A. Hippensteele 
Workshop Organizer and Editor 
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1996 COOLANT FLOW MANAGEMENT WORKSHOP 

Joe Gladden 

NASA Lewis Research Center 

Welcome to the first annual Coolant Flow Management Workshop. The purpose 
of the workshop is to review the status of work accomplished to date on a number of 
research efforts and to assess the planned work effort. The coolant flow management 
program is a response to the 1993 Turbomachinery Peer Review committee suggestion 
that LeRC pursue coupling of the internal coolant passage flow and heat transfer with the 
external gas path of the high pressure turbine through the blade/vane film cooling holes. 

This program evolved through a series of meetings with designers and heat transfer 
researchers in industry in 1994 and 1995. Modeling, turbulence models, validation data 
and grid generation were identified as topics that needed to be addressed for both internal 
passage heat transfer and film cooling. In addition, the capability to integrate the internal 
and external heat transfer solutions and optimize the results was identified as a long term 
goal. 

The workshop will cover three discreet topics: Internal Coolant Passage Heat 
Transfer, Film Cooling Heat Transfer, and Coolant Flow Integration and Optimization. 
Each session will present a combination of analytical and experimental research and will 
present both current results and future plans. There will be a group discussion at the end 
of each session for clarification of the material presented and to assess the needs for 
future research. 

NASA LeRC is developing robust multi-block schemes and multi-block merging 
techniques to improve the accuracy of the computation of heat transfer in serpentine 
passages while decreasing the time required to generate computational models. 
Turbulence models are also being assessed. Prof. Lau (Texas A & M) has been 
experimentally studying serpentine passage geometry with turbulence generating ribs at 
LeRC to provide validation data. This geometry will be modified to include bleed effects 
representing film cooling holes. Mr. Wagner (PE/UTRC) and Dr. Tse (SRA) have 
experimentally studied heat transfer in rotating serpentine passages. 

A similar effort is underway on external film cooling with the same emphasis on 
understanding the flow field around the film cooling hole and its effect on film 
effectiveness. Dr. Garg (LeRC/AYT) has developed computational techniques to 
compute the film cooling effects over the entire airfoil surface starting at the hole exit 
plane. Prof. Simon (U of Minn.) and Prof. Acharya (LSU) are studying the fundamental 
flow phenomena around the entrance to the film cooling holes to improve turbulence 
modeling. Dr. Ames (Allison) has just completed an experimental film cooling study 
using a cascade of airfoils in a high turbulence flow field. 

The ability to couple and optimize the internal coolant flow and the external film 
cooling for turbine airfoils is the goal. Research is underway at Arizona State University, 
AlliedSignal Engine Co., Ohio State University, General Electric and Penn State 
University to accomplish this goal. These programs are designed to incorporate the best 
computational tools and methodologies available and to provide validation data. 

You will be hearing about progress toward these goals during the next day and a 


half. 
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GOAL STATEMENT 


Provide an integrated enabling technology base of understanding and 
tools useful to making design decisions for the management of turbine 
coolant flow. 


References: 

"The review team strongly feels that NASA should have a leadership 
role in multi-disciplinary analysis methods. Two specific examples ... 
(ii) system that couples the internal cooling system through the film 
cooling holes to the external aerodynamics. 

♦ Turbomachinerv Peer Review Report. Sept 22. 1993 

♦ Industrial Contacts. 1994 & 1995 

♦ SAE Aerospace Atlantic Conference. 1996 


NASA / INDUSTRY IDENTIFIED NEEDS 

Internal Passage Cooling: 

♦ Modeling - wall function, NS technique 

♦ Turbulence - wall independent 

♦ Validation Data - with & w/o bleed 

♦ Grid Generation - Multi-block 

♦ Rotation Effects - 

Film-Cooling: 

♦ Methodology - code type, near-wall model 

♦ Turbulence Model - anisotropic, k-W 

♦ Validation Data - near-wall, LES, spanwise effect, high Tu 

♦ Unsteady Effects - leading edge, wall jet, cascade 

♦ Boundary Conditions - plenum effect, hole exit, conjugate problem 

♦ Grid Generation - Gridgen 

Coolant Flow Integration & Optimization: 

♦ Integration Modeling - Conjugate, NS technique 

♦ Optimization Modeling - MDO, hole/location optimization 

♦ Film-Cooling & Internal Passage Models 
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RESOURCES 


MILESTONES - Thru FY 1998/99 

Effects of Bleed Experimentally Demonstrated on Internal Passage Heat Transfer 
and Coupling Effects on Film Cooling 

Advanced Film Cooling Model Demonstrated in a Navier-Stokes Code 

Demonstrated Capability to Predict Internal Passage Heat Transfer 

Provide CFD Based Integration-Optimization Methodologies for Airfoil Coolant 
Flow 


RESOURCES 

Item FY 1994 

1995 

1996 

1997 

1998 

1999 

Total 

In-House 

$300K 

$600K 

$485K 

$445K 

$460K 

$475K 

$2,765K 

Grants 

40K 

238K 

238K 

303K 

97K 

93K 

1,009K 

Contracts 

203K 

260K 

180K 

459K 

516K 

183K 

1,801K 

Total 

$543K 

$1,098K 

$903K 

$1,207K 

$1,073K 

$751K 

$5,575K 
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Overview 

Session I: Internal Coolant Passage 

Heat Transfer 


8:30 On the Use of Structured Multi-Block Grids for Internal Coolant Flow 
Calculations David Rigby, NYMA 

8:50 Numerical Simulation of Turbine Blade Cooling Tom Shih, CMU 

9:10 Flow in Serpentine Coolant Passages with Trip Strips David Tse, SRA 

9:30 Experimental Heat Transfer in Turn Regions of Rotating Coolant 

Passages Joel Wagner, UTRC 

9:50 Experimental Multipass Heat Transfer with Trips & Bleed 

Douglas Thurman, LeRC 

10:10 Rotating Multipass Coolant Passages with Bleed Sai Lau, Texas A&M 
Break 

10:50 Group Discussion 


Heat Transfer Prediction in internal Coolant Channels 

David Rigby, NYMA, Inc. 


Objective 

• Improve heat transfer 
prediction capability in 
internal coolant passages. 

* Advance technology in this area 
to bring higher fidelity 
prediction into the design 
stage. 



A pproach 

• Multi-block approach currently being investigated. 

• Multi-block grid generation currently able to 
handle complex geometries. 

• Structured flow solver maximizes computational 
efficiency 

• Wilcox k-© turbulence model which does not 
require a reference to the distance from a wall is 
used. 


Results 



How in a channel with both ribs and bleed holes 
has been completed. 

Comparison 
of numerical 
heat transfer 
results to 
experimental 
dak is good. 


Future 

• Additional turbulence models will be investigated. 

• Progressively more complex geometries will be 
attempted. 

• Rotation effects will also be studied. 

® Grid generation and turbulence modeling remain 
the largest obstacles, followed by computational 
requirements. 
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Numerical Simulation of Turbine Blade Cooling 

Tom I-P. Shill, Mark A. Stephens, and Yu-Liang Lin 
Department of Mechanical Engineering, Carnegie Mellon University 


OBJECTIVE 

• Develop, adapt, and validate CFD codes for studying internal 
and film cooling of turbine blades. 

• Apply validated CFD codes to obtain understanding of fluid 
mechanics and heat transfer processes relevant to internal 
and film cooling and how they are affected by design and 
operating parameters. 

• Develop design strategies for more effective cooling based 
on understanding of basic flow physics. 

APPROACH 

• Internal coolant passages are currently modeled as two 
straight ducts with and without ribs connected by a tightly 
wound 1 80° bend under rotating and non-rotating conditions. 

® Film cooling is modelled by a semi-cylindrical leading edge 
with a flat afterbody with rows of compound-angle holes 
drilled through it. Computations resolve flow over leading 
edge and afterbody as well as flow in film-cooling holes and 
the internal coolant passage from which the cooling flow 
emerges. 

• Formulation is the ensemble-averaged conservation 
equations of mass, momentum (full compressible Navier- 
Stokes), and total energy closed by a “low-Reynolds 
number” k-oo model. Wall functions are not used! 

• Grid systems are multi-block structured grids - perfectly 
patched, arbitrarily patched, and overlapped. 

• Numerical method of solution utilizes flux-difference 
splitting of Roe with limiters for inviscid flux vectors, central 
for viscous vectors, and a diagonalized ADI scheme with 
multigrid for relaxation to steady state. 

• Codes used are CFL3D, MAGGIE, and RONNIE. 

• Validation is based on comparison with available 
experimental data. 


RESULTS Leading Face Trailing lace 



Ro = 0.24, Rd = 25000, Ap/p » 0.1 3 

MILESTONES 

• Showed that for a given Reynolds and Rotation number, an 
infinite number of angular speeds are possible. To fix the 
angular speed, must specify Mach number as well. 

• Showed angular speed to have significant effects on pressure 
and density variations along an internal coolant passage. 
Actual angular speeds must be simulated in order to capture 
centrifugal buoyancy correctly. 

• Showed how flow separation induced by centrifugal 
buoyancy causes the formation of additional vortical 
structures beyond the traditional two from Coriolis force. 

• Studied how secondary flows induced by Coriolis force, 
centrifugal buoyancy, 180° bend, and ribs affect each other 
and the heat transfer characteristics. 

• Studied leading-edge film cooling with injection through 
staggered rows of compound-angle discrete holes. 
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Flow in Serpentine Coolant Passages with Skewed Trips, 
David Tse, Scientific Research Assoc., Inc. 


OBJECTIVES 

• Measure, assess and analyze the flow in serpentine coolant 
passages having skewed trips. 

• Perform computational studies on the influence of Coriolis 
and buoyancy effects on the flow characteristics and heat 
transfer. 

• Analyze the velocity data in relation to the heat transfer 
measurements of Johnson, et al (1994). 

• Examine the influence of trip orientation on the flow 
characteristics and the heat transfer. 

APPROACH 

Experimental - (SRA) 

• Three components of velocity were obtained under 
stationary and rotating conditions by laser velodmetry with 
the aid of refractive-index-matching (RIM). 

Computational - (P&W) 

• The flow characteristics and heat transfer were examined 
using periodically fully-developed “conveyor-belt” boundary 
conditions. 

• Boundary conditions at the exit of the computational 
domain were used to create the inlet conditions. 


RESULTS 

• Skewed trips generate strong secondary flows and induce 
considerable streamwise development of the primary flow. 

• Rotation and skewed trips produce a swirling vortex and a 
comer recirculation zone in the straight passage. 

• Trips skewed at -45° in the outward flow passage and trips 
skewed at +45° in the inward flow passage maximize heat 
transfer. 

• The velocity measurements explain the heat transfer 
phenomena reported in Johnson, et al (1994). 

MILESTONES 

• RIM and laser velodmetry can produce highly detailed 
maps of flow field in a serpentine coolant passage. 

• RIM requires an isothermal flow field and can not include 
buoyancy forces. However, it does quantify the effect of 
the interacting Coriolis and trip driven secondary flows. 

• The overall agreement between the velodty data and CFD 
calculations is veiy good. 

• Flow, temperature and pressure fields can be assumed to be 
periodic over a streamwise repeating geometric pattern in a 
ribbed-wall passage. 
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HEAT TRANSFER IN TURNS, LET - Task 30 

Joel Wagner, United Technologies Research Center 


Objective: 

Experimentally measure, assess, and 
analyze the heat transfer within the 
internal cooling configuration of a 
serpentine coolant passage for comparison 
with CFD and numerical heat transfer 
prediction codes. 

Results: 

Thermocouple results show significant 
increases in heat transfer in the tip and root 
turn regions. Liquid crystal results show 
variation in heat transfer distribution in turns 
due to flow and thermal conditions. 

Approach: 

Conduct experiments using rotating models 
employing previously developed 
measurement techniques, including a 
transient heat transfer technique and liquid 
crystals to determine the heat transfer 
distribution throughout the model coolant 
passage. 

Milestones: 

•Design and Fabrication 
•H(0.5) Model Tests 
•H(2.0) Model Tests 
•Turn Model Tests 
•Analysis and Final Report 


Experimental Multipass Heat Transfer 
with Trips and Bleed 

Doug Thurman, Phil Poinsatte, NASA 


Objectives: To obtain simple 
\ geometry surface heat transfer 
data for internal channels with 
trips and bleed for use in 
1 validating models 

Facility: SW-6 

Approach: Transient liquid 
\crystal technique \ 

Milestones: Rib data, turn data, 
bleed data 
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LOCAL HEAT TRANSFER DISTRIBUTIONS IN 
ROTATING MULTIPASS CHANNELS WITH 

BLEED 


OBJECTIVES: 


Sai Lou, Texas A & M University 

FACILITY/RESULTS: 


Make detailed local experimental 
measurements with rotating 
multipass channel models 

• to enhance unders tandin g of effects 
of rotation, sharp turns channel 
geometry, and rib turbulators. 

• to help improve the design of these 
cooling passages. 

APPROACH: 

• Naphtalene sublimation 
experiments in rotating multipass 
channels. 


The facility consists of a rotational 
test rig, and open flow loop and 
mass transfer measurement 
equipment The results include 
normalized Nusselt and Sherwood 
number distributions. 

MILESTONES: 

• 10/99 Complete experiments with 

two-pass and three-pass channels to 
determine effect of coolant bleed. 
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Overview 

Session II: Film Cooling Heat Transfer 




1:10 Film Cooling with Anisotropic Diffusion and Hole Entry Effects: 

Experiments and Computations Terry Simon, U of Minn 

1 :30 Large Eddy Simulation and Improved Turbulence Modeling of Film 
Cooling Sumanta Acharya, LSU 

1 :50 Heat Transfer in Film-Cooled Turbine Blades Vijay Garg, AYT 

2:10 Effect of Tabs on a Jet in Cross Flow Khairul Zaman, LeRC 

2:30 Film Cooling Heat Transfer Measurements Using Liquid Crystals 

Steve Hippensteele, LeRC 
2:50 Turbine Vane Heat Transfer, Film Cooling, and Turbulence Experiments 

Forrest Ames, Allison Engine Co 

3:10 Break 

3:30 Investigation of Rotor Wake Effects on Film Cooling 

Jim Heidmann, LeRC 

4:10 Unsteady High Tu Effects on Turbine Blade Film Cooling Heat Transfer 
Using Liquid Crystals J. C. Han, Texas A & M 

4:30 Group Discussion 
5:30 Social 
6:30 Dinner 
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FILM COOLING LATERAL DIFFUSION AND HOLE ENTRY EFFECTS 

Terry Simon, University of Minnesota 


OBJECTIVES: 

To support the continued improvement of film cooling 
modeling and design by better describing the effects on film 
cooling performance of (1) anisotropic eddy transport in the 
mixing zone where freestream flow and film cooling flow 
interact and (2) the geometry of the plenum from which the 
film cooling flow enters the holes. The former will lead to 
more accurate computation of design models and the latter 
will aid in the assessment of film cooling data in the literature 
and will assist in the selection of supply plenum geometric 
features. Since detailed measurements are taken within the 
flowfield and on the film-cooled surface, the study lends a 
valuable glimpse into the physics of this complex flow. 


APPROACH: 

Experiments are conducted which simulate many features of 
the actual film cooling flow situation. In these tests, the 
freestream turbulence intensity and scale are simulated, as 
are the hole and supply plenum geometries. Measurements 
are with high-resolution traversing of fine thermocouple and 
hot-wire probes throughout the flowfield. Film cooling 
effectiveness values are found by traversing the 
thermocouple probe from the freestream down into the near- 
wall adiabatic zone to first measure the freestream 
temperature, then to deduce the adiabatic wall temperature. 
Computation is with the finite-volume numerical model with 
k-e turbulence closure modeling. 


FACILITIES: 

The facility is a flat plate with a single row of 35°, streamwise- 
odented, round film cooling holes spaced laterally at 3D and of 2.3D 
long. The delivery flow issues from an unheated combustor 
simulator section. The film cooling flow can be heated slightly for 
film cooling effectiveness measurements. 



MILESTONES: 

In HTD-VoI. 327, measurements are presented of the lateral and 
wall-normal, turbulent eddy diffusivities. These are used to 
compute the anisotropy ratio, Em,z* e m,y* I n 96-WA/HT-7, 
measurements of streamwise velocity are presented to show the 
effects on the flowfield in the interaction zone of hole length, 
freestream turbulence level, and blowing ratio. In 96-WA/HT-8, 
computational results were presented to show the effect of plenum 
geometry on film cooling performance. In a paper submitted to the 
97-IGTC, measurements are presented which show the effects on 
the flow in the interaction zone of the means by which the flow to 
the film cooling holes is channeled to the supply plenum. Presently, 
measurements are being taken with laterally-oriented holes and a test 
plate has been fabricated for compound-angle injection. Upcoming 
plans include the evaluation of performance with shaped holes and 
we are preparing to take heat transfer coefficient data. 
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TIME- AND SPACE- ACCURATE SIMULATIONS OF A FILM-COOLING JET FLOW 

Sumanta Acharya 
Louisiana State University 


OBJECTIVES 

• To perform Direct Numerical Simulations 
(DNS) and Large Eddy Simulations for a 
Film Cooling Jet How to better 
understand the large scale dynamics, their 
effects on surface heat transfer and to 
guide model development 


APPROACH 

• Higher order finite-differences for DNS 
code 

• Smagorinsky and Dynamic sub-grid scale 
stress models 

• Non-linear and other improved turbulence 
models 


RESULTS 

• Time-and space-accurate simulations reveal 

the important flow structures 



MILESTONES . 

• Codes developed for DNS, LES and RANS 
simulations 

• DNS code can handle multi-block domains 
and curvilinear grids 

• Film cooling flow and heat transfer 
simulations ongoing 
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HEAT TRANSFER IN FILM-COOLED TURBINE BLADES 

VijayGarg, AYTCorp. 


06 .4E.CIL VS 


• TO DEVELOP TOOLS & MODELS TO ACCURATELY 
PREDICT HOT GAS SIDE HEAT TRANSFER TO A 
FILM-COOLED TURBINE BLADE 


• TO MINIMIZE THE INJECTED MASS FLOW FOR A 
GIVEN COOLING EFFICIENCY 


• TO DECREASE AERODYNAMIC LOSSES INDUCED 
BY THE JET-SECONDARY FLOW INTERACTION 


RESULTS 


AERR Q A Q H 

•3-D COMPRESSIBLE NAVIER-STOKES EQUATIONS 

•VARIOUS TURBULENCE MODELS (Baldwin-Lomax 
with Mayle’s transition criterion, k-e, q-co, k-<o); 

(y+ < 1 near walls) 

•EXIT AREA OF EACH HOLE DISTRIBUTED OVER 
SEVERAL CONTROL VOLUMES 

• DIFFERENT VELOCITY & TEMPERATURE PROFILES 
SPECIFIED AT THE HOLE EXIT 

•MULTI-GRID, CELL-CENTERED FINITE VOLUME 
CODE WITH FOUR-STAGE RUNGE-KUTTA TIME 
STEPPING SCHEME, 3-D EIGENVALUE SCALING OF 
ARTIFICIAL DISSIPATION & RESIDUAL SMOOTHING 


FUTU RE- D IR ECTIO NS 


• LARGE CHANGE IN HEAT TRANSFER COEFFICIENT 
AT THE BLADE SURFACE CAUSED BY DIFFERENT 
VELOCITY & TEMPERATURE PROFILES AT THE 
HOLE EXIT 

•DIFEFRENT EFFECTS OBSERVED ON THE SUCTION 
AND PRESSURE SURFACE DEPENDING UPON THE 
BLADE & HOLE SHAPE AND HOLE LOCATION 

•BLADE ROTATION HAS A SIGNIFICANT EFFECT ON 
THE FILM COOLING EFFECTIVENESS 

• BALDWIN-LOMAX TURBULENCE MODEL SEEMS TO 
BE AS GOOD AS ANY TWO-EQUATION MODEL, AND 
IS 40-50% LESS EXPENSIVE 


• DETAILS OF IN-HOLE AND NEAR-HOLE PHYSICS 
REQUIRED UNDER CONDITIONS RELEVANT TO 
THE GAS TURBINE (DENSITY RATIO, CURVATURE 
EFFECTS, HIGH Tu, etc.) 

• EXPERIMENTAL DATA ON ROTATING FILM- 
COOLED BLADES REQUIRED 

| 

• INCORPORATE HOLE-EXIT PROFILES BASED ON 
RECENT EXPERIMENTAL DATA IN THE CODE 

• USE A MULTI-BLOCK CODE TO DEVELOP HOLE- 
EXIT PROFILES ON A ROTATING BLADE 

• INCORPORATE SIMPLIFIED FILM-COOLING MODEL 
TO SPEED-UP THE CODE 
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THE EFFECTS OF TABS ON A JET IN A CROSS-FLOW 

Khairul Zaman. NASA 


Objective: 

A tab placed on the leeward side of the 
nozzle was expected to increase jet 
penetration into the cross-flow. An 
experiment at UTRC showed 
insignificant effect. The primary 
objective of the present study was to 
confirm and explain the ineffectiveness. 

Approach: 

• Experiments in a low speed wind tunnel. 

• Conduct hot-wire measurements for 
mean velocity and stream wise vorticity 
fields. 


Results: 

• Ineffectiveness of tab, when placed on 
leeward side, is confirmed. 

• Static pressure distribution around the 
nozzle provides an explanation for the 
ineffectiveness. 

• Tab placed on the windward side of the 
nozzle is found to reduce the jet 
penetration. 

Milestones: 

• 5/97 Experiments to optimize tab 
geomery for maximum reduction in jet 
penetration. 

• 5/98 Active control for increasing jet 
penetration and mixing. 
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Film Cooling Heat Transfer w/Liquid Crystals 

2-D Film-Cooling Heat-Transfer on AlliedSignal Vane 

Steven A. Hippensteele, NASA 

Objective: Results: 

• Provide 2-D, high- * A-S design data reed 

resolution heat transfer . Scaled design points done 

map for model validation. . Cascade design started 

• Tests and Results remain 

Approach: 

• Linear vane cascade Milestones: 

• Transient Liquid Crystals . HT Data 

• Match: M, Re, c/g Den. . 1Q99-Final HT Data 

• H TCoef. and FC Effect 


Film Cooling Heat Transfer w/Liquid Crystals 
Effects of Tab Vortex Generators on Surface Heat- 
Transfer Downstream of a Jet in Crossflow 

Steven A. Hippensteele, NASA 

Objective: Results: 

• Provide 2-D, high- • Jet tests started 

resolution heat transfer • F C hole tests next 
maps for investigation. . Variation Tests and 

Results remain 

Approach: 

• Modified Sq. Duct Tunnel Milestones: 

» Transient Liquid Crystals . 2Q97-Init HT Data 

• Compare jets &FC holes; • 3Q97-Final HT Data 
w/ & w/o tabs 

• H T Coefficient Maps 
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Turbine Vane Heat Transfer, Film Cooling, and Turbulence Experiments 

Forrest Ames, Allison Engine Company 


OBJECTIVES 

o Determine the influence of high inlet turbulence on vane 
film cooling 

o Document the influence of vane film cooling on local 
heat transfer rates in the presence of high inlet 
turbulence 


APPROACH 

o Four vane ambient cascade rig 

o Mock combustor turbulence generator 

o Thin foil heat transfer technique with finite element 
analysis 

o Adiabatic effectiveness corrected for conduction 

o Hot wire boundary layer, turbulence and spectra 
measurements 

o Total pressure exit loss measurements with film 
cooling 


RESULTS 

o Suction surface film cooling results correlated well with 
flat plate data but the influence of turbulence gradually 
reduced film cooling protection 

o Suction surface heat transfer was increased in the near 
hole region. Two rows had higher augmentation than 
one 

o Pressure surface film cooling protection was strongly 
dissipated by the high turbulence 

o Heat transfer was substantially enhanced by the film 
cooling and increased with blowing ratio and with the 
number of rows 

o Single wire velocity profiles showed good spanwise 
uniformity by an X/d of 10 for the high turbulence case 
for both the single and double row cases. 


MILESTONES 

o All milestones for the contract have been completed and 
the final report has been published 
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Investigation of Rotor Wake Effects on Film Cooling 

Jim Heidmann, NASA 


OBJECTIVES : 

- use CFD to understand 3D flow physics 

- measure local film effectiveness and Nusselt number values 

- resolve chordwise, spanwise, and temporal variations 

- vary blowing ratio, density ratio, and Strouhal number 

- develop model accounting for wake-passing effects 

APPROACH : 

- modify 3D Navier-Stokes solver for unsteady film cooling 

- perform computations for experimental conditions 

- use film-cooled blade in annular turbine cascade 

- model wake passing with rotating row of cylindrical rods 

- measure blade temperatures with thin-film gauges 



MILESTONES : 


- completed unsteady computations (3rd Q ’95) 

- completed test matrix (2nd Q ’96) 

- completed data reduction (3rd Q ’96) 

- completed matrix of steady computations (4th Q ’96) 

- developed wake-affected film cooling model (4th Q ’96) 
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UNSTEADY HIGH TURBULENCE EFFECT ON TURBINE BLADE FILM COOLING HEAT TRANSFER 


PERFORMANCE USING A TRANSIENT LIQUID CRYSTAL TECHNIQUE 
J. C. Han, TurbineHetti Transfer Laboratory, Texas A&M University 


OBJECTIVES 

To study detailed distributions of heat transfer coefficients and 
film effectiveness for a gas turbine blade, including film cooling 
hole region, under the effects of upstream flow conditions: 

• Free stream turbulence 

• Unsteady wakes 

• Unsteady wakes with trailing edge coolant ejections 

and the effects of film cooling injections with air and C0 2 as 
coolants: 

• Blowing ratios 

• Coolant-to-mainstream density ratios 

• Advanced film hole shapes. 

RESULTS 

•Nusselt numbers increase with an increase in free-stream 
turbulence level, and also increase with the addition of unsteady 
wakes. Adding grid-generated turbulence to unsteady wakes 
further increases Nusselt numbers on the downstream blade 
surface. 

• The trailing edge ejection is to increase Nusselt numbers on the 
front parts of both the blade suction and pressure surfaces. 

• Film injection promotes earlier laminar to turbulent boundary 
layer transition on the suction surface and also greatly enhances 
local Nusselt numbers. 

• Unsteady wakes slightly enhance Nusselt numbers but 
significantly reduce film effectiveness on a film-cooled blade 
compared with a film-cooled blade without unsteady wakes. 

• An increase in blowing ratio increases Nusselt numbers for both 
coolants. Film effectiveness is highest at a blowing ratio of 0.8 
for air injection and at a blowing ratio of 1 .2 for C0 2 injection. 

APPROACH 

A transient liquid crystal technique has been developed for this 
study. In the experiments, the liquid crystal coated blade is heated 
to a uniform temperature, and then suddenly exposed to a pre-set j 
mainstream flow condition. A Color Frame Grabber system, 
through the RGB cameras, is used to record the time of the color 
change of the liquid crystal on the blade surface. The color change 
times are related to heat transfer coefficients and film 
effectiveness based on the solution of a 1-D semi-infinite transient 
conduction model with a convective boundary conditions. 

MILESTONE 

The effects of upstream flow conditions and film cooling injections 
have been studied. The tasks for the coming year are: 

•To design and fabricate a test blade with several rows of shaped 
film holes. These advanced film hole shape designs will be 
decided after consulting with GE for the optimum shapes. 

•To measure detailed heat transfer coefficients and film 
effectiveness for the blade with shaped film cooling holes under 
the combined effects of free-stream turbulence, unsteady wake 
and unsteady wake with coolant ejection. 
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Overview 

Session IE: Coolant Flow 
Integration & Optimization 


8:30 Multi-Disciplinary Optimization of Cooled Turbine Design Using 3-D 
CFD Analysis Yong Kim, AlliedSignal 

8:50 CFD/Expeiiment to Optimize Film-Cooling Hole Shape Yong Kim 
9:10 Film Cooling Flow and Heat Transfer Design Robert Bergholz, GE 
9:30 An Initial Multi-Domain Modeling of an Actively Cooled Structure 

Erlendur Steinthorsson, ICOMP 

9:50 Aero-Thermo-Structural Optimization 

Joe Gladden, NASA LeRC for George Dulikravich, Penn State Univ 
10:10 Break 
10:30 Group Discussion 
11:30 Adjourn 


MULTI-DISCIPLINARY OPTIMIZATION OF 
COOLED TURBINE DESIGN 

Aditi Chattopadhyay, Arizona State University 
Presented by Yong Kim, AlliedSignal Aerospace 


Objective: 

• Minimize turbine blade 
temperature and optimize coolant 
flow using Multiobjective, 
Multidisciplinary Optimization 
Procedure. 

Approach: 

• CFD flow field analysis 

• FE heat transfer analysis (blade) 

• Multiobjective, Multidisciplinary 
Optimization procedure 

- Kreis selmeier-S teinhaus er function 

- Multilevel format 


Design Variables: 

• Blade Chord 

• Surface geometry 

• Film cooling parameters 

-Hole position 

- Shape/size of hole 

- Injection angle 

-Injection rate 

Milestones: 

• FY98 Develop MDO procedure 

• FY98 2-D turbine design 

• FY99 3-D aero opt (inviscid) 

• FY99 Extend to viscous effects 
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CFD/EXPERIMENTS TO OPTIMIZE FILM-COOLING 

HOLE SHAPE DESIGN TOOLS 

Yong Kim, AlliedSignal Aerospace 


OBJECTIVE 

Develop efficient hole shape design method 
that will help achieve the contract goal of 
increasing TRIT of regional aircraft engines 
by 150 to 200°F without counterproductive 
increases in cooling flow 

RESULTS 

• CFD can be used as an effective design tool 

• Promising new hole shape was identified 

APPROACH 

MILESTONES 

• CFD as main design tool 

• CFD hole shape design efforts 

• Code calibration against experiment 

- 95% complete 

• Develop new hole shape 

• 2D cascade tests 

• Hole manufacturing trial 

- Starts in 1 51 quarter of ‘97 

• Validation test in 2D cascade rig 

• Hole shape and CFD design tool validation 
- Starts in 2 quarter of ‘98 
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Turbine Airfoil Film Cooling - Design Integration 

Robert F , Bergholz 


Objectives 

• Maximize coolant heat capacity utilization. 

• Focus on overall heat flux management via 
film and thermal barrier technologies to control 
structural thermal gradients. 

• Optimize overall film effectiveness distribution. 

• Customize film hole geometries for improved 
effectiveness with lower blowing ratios. 


Results 

• GEAE Single-Passage Cascade 

- Gas concentration / mass-transfer method. 

- T|f results for several airfoil geometries. 

• GE CR&D Transient Cascade 

- Both r| f and h measurement capability. 

- Data acquired for full-scale vane geometries. 

• GE CR&D Water T unnel 

- PLIF measurement technique. 

- Full-coverage, discrete-hole film cooling 
applicable to combustor liners. 




Approach 


Milestones 

• Acquire film effectiveness and external heat 


• 1996 

transfer coefficient data on promising film 


- Single-Passage Cascade 

cooling configurations. 


• Engine airfoil support tests. 

• Apply optimization and computational DOE 


- Transient Cascade 

methods to quantify and rank design payoffs. 


• GEAE / CR&D IR&D film tests. 

• Use unstructured grid CFD tools to predict film 


• 1997-99 

hole and external flowfields, local heat transfer 


- Single-Passage Cascade 

effects, and customized film hole geometries. 


• Engine blade / vane, showerhead, ... 
- Transient Cascade 



• AST film tests. 
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An Initial Multi-Domain Modeling of an Actively 
Cooled Structure 

Erlendur Steinthorsson, ICOMP 


Objective 

• Develop a capability for detailed simulations of 
flow and heat transfer in turbomachinery. 

• Prepare a platform for testing 
and evaluating 
turbulence models, 
film-cooling flow 
models and other ■ 
possible models 
for turbine 
cooling flows. 


A pproach 

• Multi-block grid systems for complex geometries: 

- Globally unstructured for flexibility 

- Locally structured for efficiency and accuracy. 

• Finite volume discretization of the governing 
equations. 

• Multigrid convergence acceleration for high 
computational efficiency on fine grids. 

• Algebraic and two-equation models of turbulence. 


A pplications: 

• Various internal heat transfer cases, including a 
branched duct geometry which models features 
found in coolant passages of turbine blades. 

• Flow and 
heat 

transfer in 
turbine 
cascades, 
including 
heat 

transfer on 
blade tips. 

Ongoing and future research: 

• Implement additional turbulence models. 

• Parallel computing capability. 

• Improved convergence acceleration schemes and 
preconditioning methods. 

• Automatic solution adaptive refinement capability. 




28 


AERO-TMERMO-STRUCTURAL OPTIMIZATION 
Presented by Joe Gladden, NASA LeRC for George Dulikravich, Penn State University 

Global Research Tasks Present Status 

1) Optimization of blade hot surface pressure distribution We have already developed and tested the following software: 

2) Optimization of blade hot surface shape, thermal boundary 1) Inverse design of single coolant passage shapes in 3-D coated 


conditions and inlet temperature 

3) Minimization of the number of coolant passages inside the 

blade 

4) Optimization of thicknesses of blade walls and interior struts 

5) Stress minimization in the blade material 

6) Retrieval of cold and unloaded 3-D blade shape 

Deliverables 

The tangible deliverables resulting from various phases of the 
proposed effort will include computer codes with users' manuals, 
video-tapes of the 3-D color computer graphics, semi-annual 
reports and technical papers. Computer codes to be delivered are: 

1. 3-D parametric blade outer geometry generator based on 

analytical surface patches and surface Bezier polynomials. 

2. A constrained hybrid optimizer incorporating a genetic 

evolution algorithm, a gradient search algorithm, aNelder- 
Mead sequential simplex with Rosen's projection, and our 
method for escaping from local minimas. 

3. 3-D Blade outer surface shape inverse design based on 

Navier-Stokes equations and direct transpiration concept 

4. 3-D heat conduction FEM analysis code allowing for locally 

temperature-dependant heat conduction coefficients as 
well as thermally coated and non-coated surfaces. 

5. 3-D heat conduction inverse boundary condition code for 

determining temperatures, heat fluxes and convective heat 
transfer coefficients on inaccessible surfaces. 

6. 3-D elastic deformation analysis code for coated and non- 

coated blades. 


& uncoated turbine blades, 

2) Inverse determination of temperatures and heat fluxes on 3-D 

coolant passage walls, 

3) 3-D finite element thermo-elasticity dynamic analysis of 

unsteady stresses and deformations, 

4) A hybrid constrained optimizer incorporating gradient search, 

genetic and simulated annealing. 



Axial displacements(u) in meters. 





Session I 


Internal Coolant Passage Presentations 




ON THE USE OF STRUCTURED MULTI-BLOCK GRIDS FOR INTERNAL COOLANT 

FLOW CALCULATIONS 


David L. Rigby 
NYMA, Inc. 

Brook Park, Ohio 

A.A. Ameri 
AYT Corporation 
Brook Park, Ohio 

and 

E. Steinthorsson 
NASA Lewis Research Center 
Cleveland, Ohio 


Outline 

• Objective 

• Requirements 

• Approach 

• Multi-block Reynolds averaged Navier-Stokes 

• Method of Weakest Descent used to merge 
blocks 

• Results 

• Conclusions 
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Objective 

• Production of more efficient and reliable gas 
turbines 

• Reduced design time 

• Advance technology used in coolant passage 
design 

• Higher fidelity prediction in a reasonable time 


Typical Cooled Turbine Blade 



Requirements 

• Grid generation 

• needs to be much more transparent to the user 

• Turbulence Modeling 

• flow in the coolant passage is very three 
dimensional and highly turbulent 

• using distance to a wall as a reference is 
ambiguous, and tedious to implement 

• Maximize numerical efficiency 


A pproach 


Solver 

• Reynolds Averaged Navier-Stokes 

• Multi-block approach currently being investigated 

• Structured flow solver maximizes computational 
efficiency 

• Wilcox k-co turbulence model which does not 
require a reference to the distance from a wall 
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A pproach 


Grid Generation 

• Multi-block grid generation currently able to 
handle complex geometries (GridPro) 

• Block merging done using the recently developed 
Method of Weakest Descent 


Method of Weakest Descent 

• Algorithm to merge multi-block grid system to 
minimum number of blocks 

• Increases numerical efficiency and reduces pre- 
post-processing load 

• Can also improve load balancing for parallel 
computations 
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Grade shown near candidate. Invalid patches have x. Circle 
shows candidate which is about to be removed. 


Branch Duct Topology. 283 blocks (thin lines). 
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Film Cooling Hole topology, 51 blocks 
Minimum number found was 12 



Ribbed channel with bleed holes. 283 blocks 
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*Every other grid point 


Surface Heat Transfer 






Conclusions 


• Reasonable heat transfer prediction can be 
achieved in complex geometries 

• Multi-block grid allows efficient placement of grid 
points, and efficient use of computer resources 

• Wilcox k-co turbulence model predicts heat 
transfer well, and has good numerical behavior 


Future work 

• Further development of grid generation process 
required 

• Investigation of other turbulence models 

• Pre- and post-processing still somewhat tedious 
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NUMERICAL SIMULATION OF TURBINE BLADE COOLING 


Tom I-P.Shih, Mark A. Stephens, and Yu-Liang Lin 
Carnegie Mellon University 
Pittsburgh, Pennsylvania 


Outline of Talk 

• Objectives 

• Approach 

• Results 

• Summary 


Objectives 


• Develop, adapt, and validate CFD codes 
for studying internal and film cooling of 
turbine blades. 

• Apply CFD codes to study flow physics 
associated with internal and film cooling. 

• Develop design strategies for more 
effective cooling. 
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Approach 


Formulation 

• compressible Navier-Stokes 

• low Reynolds number k-oo 
(no wall functions!) 

Algorithm 

• flux-difference splitting with limiters 

• diagonalized ADI with multigrid 

• multiblock structured grids 

(patched perfectly, patched arbitraily, or overlapped) 

Codes 

• modified CFL3D 

(rotation, initial condition, boundary conditions, ...) 

• RONNIE (patched grids) 

MAGGIE (overlapped grids) 

Validation 

Based on comparison with experimental data. 

• For internal cooling: 

Chyu (non-rotating duct with ribs) 

Wagner, et al. (rotating smooth duct) 

• For film cooling: data from Dave Bogard 
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Results 


• Internal Cooling 

• Leading-Edge Film Cooling 

Internal Cooling 


D„ = 0.5- 

>VD,= (R,+ R,)/2 = 49 
R r /D h = 41 .85 
R t / D h = 56.15 
R 0 /D h = 1.44 
Rj / D h = 0.22 
L/ D h = 14.3 
p/e = 5 



Leading Face 
Trailing Face 


^9 
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Grid System 



Rule-Based: 

• For k-co model, 5 points between y+ = 0 and y+ = 5. 

• Constant grid spacings above the wall for the first 3 
to 5 points. 

• Grid lines should be aligned with the flow direction 
as much as possible. 

• Cluster grid points to ensure that separation and 
reattachment points as well as other relevant flow 
physics are resolved accurately. 

• Grid aspect ratio can be large in boundary layers 
where grid lines are nearly aligned with the flow, but 
should be near unity in regions with the flow 
recirculates. 
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Bend 
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1-D Model Predictions of p/p jn for 
Re = 25000, Ro = 0.24, T in = 300.0 K, R m /D h = 49 
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Velocity vectors and density contours in a plane midway between the leading and trailing faces near the duct entrance 
for the “up-leg” and near duct exit for the “down-leg” for two rotation numbers. The Reynolds number for both cases is 
25 , 000 . 


Ro = 0.24 


Leading Edge 
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Trailing Edge 


Trailing Edge 
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Trailing Edge 
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Laterally Averaged Effectiveness vs. X/D 

10 F 

Surface Effectiveness at X/D = 1 .24 vs. Z/D 

— A — Computational data 

w Computational data 
— □ — Experimental data 

n o L 

— □ — Experimental data 


1 6 

8 1 

X/D 



Surface Effectiveness at X/D = 4.86 vs. 2/D 




Computational data 
Experimental data 


Surface Effectiveness at X/D = 9.98 vs. 2/D 


Computational data 
Experimental data 
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Leading Leading 



Trailing Trailing 


Leading Leading Leading 



Trailing Trailing Trailing 


Leading Leading Leading 



Trailing Trailing Tr ailing 


Velocity vectors in the cross stream plane at several strcamwise locations along the U-duct. Stations 1, 2, 3, and 4 are 
respectively 1.45D h , 3.65D h , 5.15D h , and 7.4D h from the duct entrance. Stations 10, 1 1, and 12 are respectively 
3.45D h , 6.45D„, and 14.3D h measured from the end of the bend (i.e.. Station 9). Ro = 0.24, Re = 25,000. 
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Convex 
















Results 


• Internal Cooling 

• Leading-Edge Film Cooling 

Film Cooling 



IEA01NS EOGE - FflOWT VIEW 0E7AH. 



10 hoi e* on better# row 



u„ = 10 



m/s, = 1 atm, T a = 275° C, I„ = 05% 


(p U )c 

(pu). 


= 2.0, m c = 145 g/s 
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Summary 


Internal Cooling: 

• For smooth U-duct, specifying Re & Ro is 
inadequate. The actual rotational speed 
must be simulated. 

• For U-duct with ribs, more study needed on 
the effects of rotational speed. 

• Results showed secondary flows induced by 
Coriolis, centrifugal buoyancy, inclined ribs, 
and 180 ° bend, and how they interact and 
affect heat transfer. 

Film Cooling: 

• Need low-Reynolds no. models that can 
account for shear layers at the wall and at 
jet/freestream interaction region. 

• Also need turbulence models to account for 
anisotropic spreading (lateral versus 
normal). 




FLOW IN SERPENTINE COOLANT PASSAGE WITH TRIP STRIPS 


David Tse 

Scientific Research Associate, Inc. 
Glastonbury, Connecticut 


BACKGROUND AND INTRODUCTION 


• In advanced gas turbine engines increased temperatures, stage pressure 
ratios and rotor speeds are used to increase thrust/weight ratio and to 
reduce specific fuel consumption. 

• Efficient internal cooling is essential to maintain structural integrity. 

• Rotation gives rise to Coriolis and buoyancy forces which can 
significantly alter the local heat transfer in the coolant passage. 

• A better understanding of interaction of Coriolis, buoyancy and trip 
induced secondary flows and the capability to predict heat transfer 
response to these effects is necessary for achieving efficient cooling. 


• The complex coupling of Coriolis and buoyancy forces - many 
investigators to study the heat transfer characteristics of rotating cooling 
passage. 

• Hajek, et al. (1991), Wagner, et al. (1991), Yang, et al. (1992), Mochizuki, 
et al. (1994) - experimental studies on heat transfer performance in 
rotating serpentine passages of square cross section with smooth walls. 

• Significant increase in heat transfer at the turns and there were 
considerable differences between inward and outward flow in the straight 
passage. 

• Velocity field - plays an important role in convective heat transfer. 
Velocity characteristics for flow in a passage smooth wall: Tse, et al. 
(1994) - NASA CR-4584, Tse and McGrath (1995), McGrath and Tse 
(1995). These results explain many of the heat transfer phenomena noted 
by previous investigators. 
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BACKGROUND AND INTRODUCTION 


• In practice, cooling passages contain trips, which create secondary flow to 
augment heat transfer. 

• Taslim, et al. (1991a, 1991b), Wagner, et al. (1992), Prakash and Zerkle 
(1993), Johnson, et al. (1993, 1994) indicate that the secondary flow 
produced trips and trip orientations have large effects on heat transfer. 

• This presentation deals with flow in a serpentine passage with ribbed 
walls. 

• To explain the results of Johnson, et al. (1994) 

1. The large decrease in heat transfer in the inlet region, measured in 
models with smooth walls, does not occur in models with trips. 

2. The effects of density ratio and flow direction for models with trips 
are less than for models with smooth walls. 


OBJECTIVE 


• Benchmark quality velocity measurements were acquired by LDV at 
Re = 25,000 and Ro = 0.24 in a rotating serpentine passage having skewed 
trips. 

• To assess the influence of Coriolis and trip induced secondary flows on 
the velocity characteristics. 

• To explain the heat transfer phenomena observed in Johnson, et al. 
(1994) -ASME J. Turbomachinery, Vol 116. 

• For validation of computational procedures. 

• Results presented here are reported in detail: Tse and Steuber (1996) - 
NASA CR-xxxx (In preparation); Tse and Steuber (1996a); Steuber, et al. 
(1996). 
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APPROACH - 1: Experimental (SRA) 


Velocity measurements were obtained with the aid of refractive-index- 
matching (RIM). 

RIM - advantage 

• Avoids defraction of laser beams at the solid liquid interface 

• Allows measurements in complex geometries with multiple obstructions 
and curved internal surfaces 

• Density and viscosity are close to those of water 

• Can match Re and Ro of most turbomachinery applications 

RIM - limitation 

• RIM fluid is a liquid 

• Cannot address buoyancy 


APPROACH - 2: Computational (Gary Steuber of P&W) 

• Flow and heat transfer characteristics were examined using periodic full- 
developed “conveyor-belt” boundary conditions. 

• Boundary conditions at the exit of the computational domain were used to 
create the inlet conditions. 
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SRA Rig View from outside of model ’’looking in” 



Measuring Locations of Passage 1. 
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RESULTS (SUMMARY & COMPARISON WITH SMOOTH WALL) 
Three components of velocity - Re = 25,000 and Ro = 0.24 

• First passage (radially outward flow passage) 

Ribbed wall Smooth wall 

Entrance effect - 1.0 D Entrance effect - 6.5 D 

Coriolis & trip secondary flow Coriolis effect 

• First turn (Ribbed wall & Smooth wall) 

Effect of the turn - Swirl 

• Second passage(radially inward flow passage) 

Ribbed wall Smooth wall 

Effect of the turn - 1.0 D Effect of the turn- 4.0 D 

Coriolis & trip secondary flow Coriolis effect 

• Second turn (Ribbed wall) 

Effect of the turn - Swirl (weaker than the first turn) 


y/H « 1 



1 


y/H = l 



Velocity vectors in the first passage. 
1 yJD = 1 2 x/D = 7 
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ID downstream of T1 ID upstream of T2 


Streamwise velocity contour 
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V-Component and W-Component Velocity Contours - Stationary Channel Location A2 
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(SR A Data) W-Velocity Component 
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V-Component and W-Component Velocity Contours 


Rotating Channel Location D3 




(a) Compressible Flow Channel Simulation (b) Compressible Row "Conveyor- belt” Simulation 

Rotating Channel (at i- Plane = 760) Rotating Channel Segment 


Comparison of Secondary Flow Velocity Vector Fields — Channel vs “Conveyor-Belt” Simulation 


CONCLUSIONS - 1: Stationary measurements (Re = 25,000 and Ro = 0.0) 

• Skewed trips generate strong secondary flows. Outward flow and trips 
skewed at -45° lead to a two-vortex structure with upward velocity of 
0.6 U b in the center and a downward velocity of -0.95 U b at 0.09 D from 
the wall. When both trip orientation and flow direction are reversed, i.e., 
inward flow and trips skewed at +45°, the vortices circulate in the same 
direction. When either trip orientation or flow direction alone is 
reversed, the vortices circulate in the opposite direction. 

• The streamwise development of primary flow is considerable, as a 
consequence of the strong secondary flow. Skewed trips transform a plug 
flow to two regions of high velocity, velocity exceeding 1.50 U b close to the 
side walls within 12 D of entrance. 
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CONCLUSIONS - 2: Rotating measurements (Re = 25,000 and Ro = 0.24) 

• With clockwise rotation, outward flow and trips skewed at -45° generate 
counter-clockwise swirl and a corner recirculation zone at the upper 
corner of the leading side. Inward flow and trips skewed at +45° generate 
clockwise swirl and a corner recirculation zone at the upper corner of the 
trailing side. Outward flow and trips skewed at +45° generates clockwise 
swirl and a corner recirculation zone at the lower corner of the trailing 
side. 

• With trips skewed at +45°, the secondary flows remain qualitatively 
unaltered as the cross-flow proceeds from the first passage to the turn. 
The flow characteristics at these locations differ when trips are skewed at 
-45°. Changes in the flow structure are expected to augment heat 
transfer, in agreement with the heat transfer measurements of Johnson, 
et al. (1994). Trips skewed at -45° in the outward flow passage and trips 
skewed at +45° in the inward flow passage maximize heat transfer. 


CONCLUSIONS - 3: Computational (Re = 25,000 and Ro = 0.24) 

• The overall agreement between the SRA velocity data and the P&W CFD 
calculations is very good. 

• The agreement between the “conveyor-belt” segment and the full skewed 
trip strip channel simulation is excellent. 

• Flow, temperature and pressure fields in a ribbed-wall passage can be 
modeled by the assumption of a periodic repeating pattern in the 
streamwise direction. 

• The use of “conveyor-belt” segment reduces the computational domain 
from ~1, 700, 000 grid points to -100,000 grid points (less than 6% of the 
requirement for full channel simulation). 
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MILESTONES 


Experimental 

• RIM and laser velocimetry can produce highly detailed maps of flow field 
in a serpentine coolant passage. 

• RIM requires an isothermal flow field and can not include buoyancy 
forces. However, it does quantify the effect of the interacting Coriolis 
and trip driven secondary flows and provides one more piece to the 
puzzle. 

Computational 

• Flow, temperature and pressure fields in a ribbed-wall passage can be 
modeled by the assumption of a periodic repeating pattern in the 
streamwise direction. 

• The use of “conveyor-belt” segment reduces the computational domain 
from -1,700,000 grid points to -100,000 grid points (less than 6% of the 
requirement for full channel simulation). 


FUTURE WORK 


• SRA velocity data show that the flow in the first passage is characterized 
by a swirling vortex and a corner recirculation zone but the P&W CFD 
calculations show a third recirculation bubble in the inter-rib region. 

• The largest difference between the experimental data and the numerical 
simulation is in the inter-rib region. The differences could be in part due 
to the sparse nature of the data set near the wall. The resolution of the 
existing LDV system in the inter-rib region is inadequate. 

• Improved resolution can be achieved by microscopic imaging and PIV. 
PIV utilizes a laser sheet and resolves the diffraction problem associated 
with the crossing of two laser beams. 

• Microscopic imaging and PIV has the potential to investigate air flow and 
address experimentally the effect of buoyancy. 
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ASSESSMENT OF HEAT TRANSFER COEFFICIENTS IN TURNS OF TURBINE 

BLADE COOLANT PASSAGES 


Joel Wagner 

United Technologies Research Center 
East Hartford, Connecticut 


Objective 


• Experimentally measure, assess, and 
analyze the heat transfer within the 
internal cooling configuration of a 
serpentine coolant passage for comparison 
with CFD and numerical heat transfer 
prediction codes. 


Approach 


• Conduct experiments using rotating models 
employing previously developed 
measurement techniques, including a 
transient heat transfer technique and liquid 
crystals to determine the heat transfer 
distribution throughout the model coolant 
passage. 
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Results 


• Thermocouple results show significant 
increases in heat transfer in the tip and root 
turn regions. Liquid crystal results show 
variation in heat transfer distribution in turns 
due to flow and thermal conditions. 
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Douglas Thurman and Philip Poinsatte 
NASA Lewis Research Center 
Cleveland, Ohio 


Testing is currently in progress to acquire simple geometry surface heat transfer 
data for internal channels with trips and bleed holes which can be used in the 
development and validation of models. The transient liquid crystal technique is 
used on a simple multipass model with rectangular channels and normal ribs. 
Normal bleed holes are located on the floor of the model in the first channel. 
Each hole is attached to a flow meter, allowing various bleed flow rates to 
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HEAT TRANSFER LOCAL DISTRIBUTIONS IN ROTATING MULTIPASS 

CHANNELS WITH BLEED 


S.C. Lau, C.W. Park, R.T. Kukreja, and M. Kandis 
Texas A&M University 
College Station, Texas 

OBJECTIVES: 


The main objective is to make available 
detailed local experimental data with rotating 
multipass channel models: 

♦ to enable better understanding of the effects of 
rotation, sharp turns, channel geometry, and 
rib turbulators on the heat transfer 
distributions in cooling serpentine passages in 
modem gas turbine blades; and 

♦ to help improve the design of these cooling 
passages. 


HEAT TRANSFER IN SERPENTINE 
PASSAGES IN TURBINE BLADES: 

♦ Row Parameters — Reynolds number. Rotation number, 
and buoyancy parameter relate viscous, inertia, Coriolis, and 
buoyancy forces in the flow. 

♦ Thermal Parameters — Heat transfer distribution also 
depends on Prandti number and thermal boundary 
conditions. 

♦ Geometric Parameters: entrance geometry, passage cross 
section, variation of cross section along the main flow 
direction, wall orientation with respect to rotation axis, turn 
configuration, and inward or outward main flow direction. 

♦ Turbulator Parameters: rib cross sectional shape, size, and 
pitch; rib angle and orientation with respect to main flow 
direction; and rib configurations on opposite walls. 
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SELECTION OF PRIOR PUBLICATIONS: 


Turbulent flow in unheated channels : 

Johnston et al. (1972), Wagner and Velkoff 
(1972), Rothe and Johnston (1979), and Tse and 
McGrath (1995) 

Turbulent heat transfer in straight channels : 

Mori et al. (1971), Morris and Ayhan (1979), 
Clifford et al. (1984), Harasgama and Morris 
(1988), Soong et al. (1991), Taslim et al. (1991a, 
1991b), El-Husayni et al. (1994), Hwang and 
Kuo (1994), and Kuo and Hwang (1996) 


SELECTION OF PRIOR PUBLICATIONS: 


Turbulent heat transfer in multipass channels : 

Hajek et al. (1991), Wagner et al. (1991a, 1991b, 
and 1992), Johnson et al. (1993, 1994), Yang et 
al. (1992), Mochizuki et al. (1994), Han and 
Zhang (1992), Han et al. (1993), Parsons et al. 
(1994), Zhang et al. (1995) 

Numerical Studies : 

Prakash and Zerkle (1992), Dutta et al. (1994a), 
Bo et al. (1995), Tse and Steuber (1996), and 
Shih et al. (1996) 
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APPROACH: 


Naphthalene sublimation experiments are 
conducted to obtain the local mass transfer 
distributions in rotating multipass channels. 
By applying the heat/mass transfer analogy, 
the distributions of the normalized Nusselt 
number may be determined from the 
measured distributions of the normalized 
Sherwood number. 


NAPHTHALENE SUBLIMATION 
TECHNIQUE: (Goldstein and Cho, 1995) 

♦ Same governing equations for heat and mass 
transfer — negligible viscous dissipation and 
compression work in energy equation; no 
coupling between temperature and 
concentration fields; negligible property 
variation; same Pr t and Sc t (Simpson and 
Field, 1972). 

♦ Boundary condition — naphthalene surface is 
analogous to isothermal surface. 
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NAPHTHALENE SUBLIMATION 
TECHNIQUE: (Goldstein and Cho, 1995) 

♦ Negligible disturbance of velocity field by 
sublimation — small change of surface shape 
and low normal velocity at naphthalene 
surface. 

♦ Negligible surface temperature change due to 
latent heat of sublimation. 

♦ Heat/ mass analogy: Nu/Sh = (Pr/Sc) n . If Nu r ^ 
cc Pi* and Sh re f qc Sc 11 , 

Nu/Nu re} = StySh ref . 


FACILITY/RESULTS: 

♦ The experimental facility consists of a 

rotation test rig, an open air flow loop, and 
mass transfer measurement equipment, in an 
air-conditioned laboratory. Typical results 
include the local distributions of Nu/Nu.f d or 
Sb/Shf d on the leading and trailing walls of a 
test channel and the distributions of the 
spanwise or regionally averaged Nu/Niif d or 
Sb/Shf d along the main flow direction on the 
leading and trailing walls. 
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MILESTONES: 


Results have been obtained for rotating 
two-pass square channels with smooth walls 
and with 90° and 60° ribs of various sizes on 
the leading and trailing walls: Reynolds 

number ranged from 5,500 to 14,500, and 
rotation number was varied up to 0.24. 
Attention was focused on the sole effect of 
rotational Coriolis forces under zero buoyancy 
and isothermal wall conditions, and its 
coupling with the effects of turn-induced and 
rib-induced secondary flows. 


results: 


♦ Re and Ro effects for a smooth channel and a 

O 

channel with 90° ribs 

♦ blockage effect of 90° ribs on leading wall on 
trailing wall heat transfer 

♦ blockage effect of 90° ribs on trailing wall on 
leading wall heat transfer 

♦ rib angle effect with 90° and 60° ribs 

♦ rotation-turn coupling effect and rotation-rib 
coupling effect on local heat transfer 
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RESEARCH PLAN: 


Experiments will be conducted with two- 
pass channels and three-pass channels to 
study the effects of coolant extraction through 
bleed holes, channel geometry, and channel 
orientation (with respect to the rotation axis) 
on the local heat transfer distributions on the 
leading and trailing walls of the channels. 
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GROUP DISCUSSION - DECEMBER 12, 1996 

SESSION I. INTERNAL COOLANT PASSAGE PRESENTATIONS 


Kestutis C. Civinskas - Facilitator Phil Ligrani - Scribe 

Clarifications and questions. 

Status of work. 

Assessment of what needs to be done. 

Rob Norton 

How many grid points - full internal cooling geometry? How many points? 
Dave Rigby [Answer] 

Millions. 

600,000 for results shown today. 

Depends on turns. 2.4 million to do a 180° turn. 

Inside & outside blade, 5 million. 

Hopefully with models, you will not have to grid every feature. 

Structured grid takes many manhours to generate grid. 

How about unstructured methodology to get a grid? 

Why have fast solvers when so much time goes into grid generation? 
[Answer] 

Grid line generation, structured, unstructured - long way to go. 
Unstructured, resolve viscous layers - not automatic. 


Joe Gladden 


Need impact in the ‘98-‘99 timeframe. 

Why structured direction? 

Unstructured takes longer. 

Multi-block question - does it take long? 

Joel Wagner 

Unstructured - where do we need to go. Number of grid points drop a lot for adaptive 
unstructured meshes. 
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Tom Shih 


Grid generation - not just geometry, but also physics. 
Boris Glezer 


How to model? 

Two processes: 

Academic - detailed flow predictions. 

Industry - reasonable/short term. 

More years needed to get details. 

Rotation is a different story - need simplified correlations. 
Rotation is first priority for industry. 

Bob Simoneau 


What level of spatial resolution (for heat transfer) is important? 

People doing research, what is really needed to know is: what are best estimates of what is 
needed? 

Kinds of resolution and locations of detail needed. 

More CFD, unstructured direction. 

Bob Bereholtz 

Simoneau’ s point is important. 

Probabilistic design systems - what is required precision for a particular type of computation. 
Given natural uncertainty: Wagner picture of a turbine blade with + values is helpful. 
Example: combustor profile is an important issue regardless of grid points. 

Need balance precision. 

Boris Glezer 

What level of resolution? 

Spacing? 

Internal channel different from outside structure? 

Third order of interest - very discrete grid - predicted results. 

Yong Kim 

Did Rigby look at any particular turbulence models other than k-omega? 
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Dave Riebv 


Not really. We did try the Baldwin-Lomax model on the branched duct and the results were 
pretty good. 

Ali Ameri 


Choose particular model for generality. 

Joel Wagner 

CFD used to calculate losses. 

Takes large number of iterations to get loss. 

Need times 10 normal grids to get losses and maybe even more to get heat transfer. 
Dave Riebv 

If you want heat transfer, it takes longer to get global flow features. 

Joel Wagner 

Accuracy required, Simoneau. 

50°F means a factor of 2 in live. 

Need surface temperatures to 5 to 10° F of actual for temperatures of 2000-3000° F. 
Bob Simoneau 


Spatial resolution also important. 

How locally universal is flow physics. 

Measured very little of local gas temperatures. 

All were referenced to some inlet conditions. 

If you have total temperature distribution, how right is local heat transfer distribution. 
Computationally, there might be some help. 

Joel Wagner 

For heat transfer coefficients, we struggled a long time because of using the bulk mean 
temperature. 

Likewise, for laminar and turbulent calculations over an airfoil to predict transition, how do you 
get momentum thickness on a 3-D airfoil flow? What is the driving potential? We struggled 
with that. 

For transient tests its even more pronounced. 
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Kaz Civinskas 


Are there questions for other speakers ? 

Metal temperatures, hot. 

Do we need 10 times the number of grid points or we will miss some physical phenomena? 
Boris Glezer 


Figure with rotation? 

Joel Wagner 

Yes, P&W blade; inside passage - correlations for rotating effects. 

Mine-King Chvu 

Why physics - fewer points. Why is temperature lower? 

Joel Wagner 

Passages have low pressure and high pressure sides that are not related to the external pressures. 
Heat transfer is up on the high pressure side of an internal passage which might be on either side 
of the blade. 

One spot / a drop / different flow directions. 

A flow divider wall is there. 

Jim VanFossen 


Liquid crystal measurements - what is best reference temperature for a gas? 

What does engine designer need? 

Joel Wagner 

Use local coolant flow, bulk mixed flow mean for current internal coolant flow design systems. 
Aerodynamic and internal cooling calculations are not separate anymore. 

Current design systems use the local bulk mean temperature. 

Copper model, hot. 

Rob Norton 


Cannot scale Nusselt number to inflow conditions, you must use local temperatures. 
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Mine-King Chvu 


Local bulk mean temperature, ways to correct now, papers in next IGTI Congress - there are 
ways to correct for this nonlinear change. 

Experiment - more convenient to use inlet temperature. 

Steven Hinnensteele 

Thurman and Poinsatte measured data recently in a large scale model having 3 pass internal 
cooling flow with ribs. They recorded temperatures as function of time and location throughout 
the transient liquid crystal test. The temperatures looked almost linear with location for many 
cases. 

J. C. Han 


Transient large crystals - two parts. Thermocouples in my experiment in different locations. 
High Reynolds numbers - 50,000/60,000: almost no change of temperature with time. 

Dave Rigby 

Spanwise integration of bulk temperature showed very linear results. 

Much numerical data - post processing of numerical data. 

Experimentally, depending on where global energy addition, can be linear. 

Curved portions different. 

Isothermal surface conditions. 

Bob Bergholtz 

What temperature to use in blade design? 

What is dovetail condition - designer, look at various changes in engine, seal conditions, etc. 
Either. Design blade or reduce uncertainty in blade life. 

Thermal and mechanical interactions are driver to these computations. 

Look more at engine in reducing design uncertainty - interesting computationally. 

Rob Norton 


No such thing as a perfect engine. 
Must design for good and bad engines. 
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Ron Bunker 


Need to reduce design uncertainty, and increase accuracy in prediction methods. 

Turbine designs - designs pushing to other limits, i.e. no film cooling at all. 

Accuracy on CFD predictions will be tighter. 

Still need a lot more detail. Goal of coolant flow management is to take out as much coolant as 
possible. 

Experimental basis - need tighter knowledge of uncertainty. 

Satish Ramadhvani 

High thermal stresses in the vicinity of the holes. 

Large uncertainties in approach flow. 

Need resolution in the computations for flow details. 

Rob Norton 

80% of the cost of the engine (or components) is decided in the first 20% of design cycle. 

Feed data into correlations which designer needs for estimates. 

Bob Simoneau 


Bleed hole question on internal channels: just one more variable in a complicated problem. 

What bleed conditions to test should come from industry. 

Numerical people - where to put their energies to address this? 

Delighted to see questions about how can these people direct research so that results are obtained 
quickly. 

Sumanta Acharva 

Must define a fine mesh, not coarse grids, to get fine flow details. 

Bob Bergholtz 

Promote - not a plea for in-precision. Detailed accuracy is needed. 

Computation - do simulation, low order internal passage configuration. Then compute blade 
flow - prescribed flow characteristics of flow and coolant: lower order resolution of internal flow 
passage. 

Very interesting to see; help to reduce uncertainties. 

David Rigbv 

From simulation, several calculations. Changes give reliable trends. 

If you don’t do high resolution calculation - not same result as low resolution calculation. 
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Bob Bereholtz 

Careful to put precision emphasis correctly. 
Tom Shih 


Geometries - internal cooling not that complicated. 
Solver/structured and unstructured meshes. 

Rob Norton 


Not interested if unstructured; less accurate if we have good test data. 
Erlendur Steinthorssen 


Grid structure, grid points are a concern. 

Need more resolution with an unstructured grid. 

Detailed simulation/more viable for structured grid. 

[To Dave Rigby] Pitch-to-height ratio covered but not correct value. Then vary all cases over a 
weekend. 

Dave Riebv 

Grids done ovemight/topology is established if doing grids of holes. 

Rob Norton 


Working on CAD systems to change geometry is separate from grid generation. 
Unstructured meshes in short term. 

Erlendur Steinthorssen 

Grid generation and block structures done simultaneously. 

Adaptive refinement can also be done in structured grids. 

Proved in a workshop some time ago. 

Structured grid codes better able to calculate some flow features. 

Rob Norton 


Better ways to estimate errors needed - rather than just depend on gradients. 
Erlendur Steinthorssen 


Reference temperatures and internal flows are important - need more than pieces of information. 
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Yong Kim 

[To David Tse] Conclusion: is heat transfer better with -45° and +45° ribs? 
David Tse 

Combination of rib angle and rotation. 

Ribs are important. Are positive in one channel and negative in another channel. 
Detailed explanation of relations between heat transfer and flow patterns. 

Kaz Civinskas 

Comments. 

Adjourn. 
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Session II 


Film Cooling Presentations 




FILM COOLING LATERAL DIFFUSION AND HOLE ENTRY EFFECTS 


Terry Simon 
University of Minnesota 
Minneapolis, Minnesota 


PRESENTATION OUTLINE: 


• INTRODUCTION: 

• OBJECTIVES: 

• EXPERIMENTAL FACILITY & 
INSTRUMENTATION: 

• MEASUREMENTS: 

• Anisotropy of Turbulent Transport 

• Velocity and Temperature Fields to Document: 

• Velocity Ratio Effects 

• Free-Stream Turbulence Effects 

• Hole Delivery Length Effects 

• Effects of Delivery Plenum Geometry 

• CONCLUSIONS: 


OBJECTIVES: 

Support film cooling modeling and design by describing the 
effects on film cooling performance of 

(1) anisotropic eddy transport, and 

(2) the geometry of the plenum from which the film 
cooling flow enters the holes. 

In the process, we documented the effects of: 

(1) Freestream Turbulence 

(2) Hole L/D 

(3) Blowing Ratio 

• Will lead to more accurate computation with design models, 

• Will aid in the assessment of film cooling data in the 
literature, and 

• Will assist in the selection of supply plenum geometric 
features. 

Also, lends a glimpse into the physics of this complex flow. 
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Nozzle 


Turbulence 

Generator 


Settling 

Chamber 



Flow From 
Secondary Fan 


Film Cooling 
Supply Plenum 



Flow From 
Secondary Fan 


High Turbulence Facility 
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Film Cooling Test Facility 



Freestream Facilities: 

• High FSTI: 12% at 10.8 m/s 

• Low FSTI: 0.5% at 10.8 m/s 


Film Cooling Configuration 


• One Row of 1 1 holes 

• 35° Streamwise Injection 

• L7D=7.0 and 2.3 

• Hole Diameter = 1 .9 cm 

Coolant/Freestream 
Parameters and Ratios 


Freestream 
Supply Row 
► 


• Velocity Ratio = 1.0 (0.5) 

• Hole Reynolds Number= 
1 3,000 (6,500) 

• Density Ratio = 1.0 


Upstream Test Downstream 
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Film Cooling Geometries 


Frees tream 



Co-Flow Delivery 


Velocity distribution at x/D= 5.0, z=zo 

(H=12% FSTI, L=0.5% FSTI) 



0 3 6 9 12 

y/e 
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, Anisotropic turbulence model 



Fig. 1 1 Cooling effectiveness distributions predicted by the isotropic 
and anisotropic turbulence models 
( from Sathyamurthy and Patankar, 1992) 


129 







Effectiveness 


COMPUTATION 



Comparison of measured and predicted values of film-cooling effectiveness 
( from Sathyamurthy and Patankar, 1992) 

M=0.5 
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on /n 


Velocity Distribution at VR=0.5, x/D=5, 
FSTI=12% for three z-planes 



y/0 


Boundary layer velocity right behind holes 
is strongly affected by film cooling flow. 


Turbulent diffusion of momentum: 


y-transport: 

z-transport 



normal to the wall 
lateral direction 


anisotropy 
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£M,z/£M,y 


em,z = - u' w' / (dU/dz), 
em,y = -uV / (dU/dy) 


Eddy diffusivity ratio at different turbulence intensity 



COMPUTATION WITH STRONG ANISOTROPY 
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CONCLUSIONS THUS FAR: 

• The eddy diffusivity ratio, > is in the range 4-15. It 

begins high and asymptotically decays, reaching about 4 by 
x/D=5.0. Anisotropy is not significantly dependent on 
freestream turbulence or blowing ratio. 

• Computational evidence suggests that use of these measured 
anisotropy values improves predictions of effectiveness 
values between the holes. 

• With low VR, the film cooling flow acts like a blockage to 
the main flow. The blockage effect is greatly reduced after 
10 hole diameters (but the vorticity imparted on the flow 
persists). 



Measurements 


Hole-Exit Profile 
Measurements 

• Document Differences in 
Effective Velocities and 
Turbulence Due to Hole- 
Length 

Flowfield Measurements 

• Document Distributions 

of “Streamwise” Velocities 
and Local Turbulence 
Levels at x/D=2.5 and 5.0 

• Compare Distributions for 
Different Cases to Quantify 
Effect 



x/D=0 x/D=2.5 x/D=5.0 
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on/n 



Mean Velocity 

Long L/D, VR=1.0, x/D=2.5, FST1=12% 
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1.6 1.4 12 1 03 06 * 04 02 ~~ 0 '”- 0.2 - 0.4 

Z/D 

% Difference in Normalized Mean Velocity 
Long vs. Short L/D 
(VR=1 .0, x/D=2.5, FSTI=1 2%) 


T"“ I i ' »' 1 * 4 I i i i 



2/D 


% Difference in Normalized Mean Velocity 
High vs. Low FSTI 

(VR=1 .0, x/D=2.5, Short L/D) 
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(a) VR=0.5 


Adiabatic Effectiveness for 
Short-Hole Injection at 
Two Velocity Rati os. 


FSTI = 12% 




Supply Plenum 


Hole Profile Measurements 
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ftole-Exit Measurements 

Normalized Mean Effective 
Velocities along Hole Centerline 



CONCLUSIONS THUS FAR: 


• Short hole injection results in "jetting" characterized by 
further penetration of the film cooling flow into the 
mainstream. It has a profound effect on the mixing zone 
between the two flows. 

• The flow between holes is more strongly influenced with 
short hole injection than with long holes. This effect 
diminishes with elevated FSTI. 
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U/Ufs 


Hole Centerline Velocity Profiles 
FSTI=12%, y/D=0.0027 


* I A ATT 

I T ♦ ♦ 


Short, VR«1 .0 1 
(^nterflowrA/R-4-A 
Cio-Flow, VR=i.O 


Span of Film Cooling Hole (x/W) 



Short-Hole VR=0.5 

Short-Hole VR=1.0 
Counter-Flow VR=0.5 

Counter-Flow VR=1 .0 
Co-Flow VR=0.5 
Co-Flow VR=1 .0 


0 


2 






Flow Streamlines for the Various Injection Geometries: 

(a) . Short hole injection with an open plenum. 

(b) . Short hole injection with counter-flow approach. 

(c) . Short hole injection with co-flow approach. 

(d) . Long hole injection. 


CONCLUSIONS FROM TfflS 
PHASE: 


• The plenum geometry has a significant effect when the 
plenum height is two diameters or smaller. Sensitivity to the 
way in which the flow is delivered to the holes rises with 
VR. 

• Computational evidence is consistent with the experimental 
findings. 

• The plenum geometry is of secondary importance relative to 
the hole length. 
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NUMERICAL SIMULATION OF A FILM COOLING JET IN CROSSFLOW 


Sumanta Acharya 
Louisiana State University 
Baton Rouge, Louisiana 

Background 

■ The hydrodynamics of a 
film-cooling jet in 
crossflow is rather 
complex. Near field is 
dominated by coherent 
vortical structures.These 
include: Jet-shear layer 
vortices, Counter-rotating 
vortex pair, Horse-shoe 
vortices and Wake 
vortices. Need to 
accurately predict these 
dynamical structures. 

■ In jet and mixing-layer flows, the near-field 
dynamical structures control the mixing and 
entrainment behavior. 

■ The dynamics of these structures can not be 
predicted by time-averaged predictive 
procedures. 




Film Cooling Jet: Measurements 

■ Measurements of Ajersch et al. (1 995) 

■ Square coolant hole$(12.5 mmm); Jet Velocity=5.5 
m/s; Blowing ratio-0.5,1 ,1 .5; Ninety-degree simple 
injection; Evidence of dynamical effects. 



From Ajersch et al (1 995) 

■ Reynolds-averaged predictions do not compare well 
with measurements 


Film Cooling: Measurements 


■ UTRC data; Circular film 
cooling hole(22.73 mm 
dia); Free stream a * 
velocity=15.25m/s; 

Blowing ratio=0.5,1, 2; 
Simple injection at 35 « 

degrees. 

■ Reynolds-Averaged 
predictions over-predict M 
jet penetration due to 
inaccurate representation 
of large scale effects. 
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Film Cooling Effectiveness 


■ Due to overprediction 
of jet penetration, film 
cooling effectiveness 
is poorly predicted. 


■ Due to overprediction 
of jet penetration, film 
cooling effectiveness 
is poorly predicted. 



FROM UTRC DATA AND REPORT 
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Goal of This Study 

■ To accurately predict the near-field large scale 
structures, and to model the smaller-scales. 
Goal is to correctly predict the jet penetration so 
that calculations of film cooling effectiveness are 
meaningful. 

■ Toward this goal, Direct Numerical Simulation 
(DNS) and Large Eddy Simulation (LES) will be 
undertaken. Comparison with RANS predictions 
will guide the evaluation and improvement of 
turbulence models 


RESEARCH PLAN 


■ Ajerschetal 

■ DNS 

■ LES 

■ RANS 

■ Improved understanding 
and predictive capability 

■ Completed DNS code 

■ Ongoing DNS simulations 
and LES simulations 


■ UTRC Data 

■ LES 

■ RANS 

■ Ongoing RANS 
simulations 

■ LES simulations to follow 
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DNS Code Development 


■ INS-3D 

■ Pseudocompressibility with 
pseudo-time iterations 

■ Fifth-order accurate 
approximation for the 
convective terms 

■ Fourth order accurate 
approximations for the 
diffusion terms 

■ Second order accurate 
temporal discretization 

■ Body fitted coordinates 

■ Fourth-order approximation 
for all metrics 


■ INS-Poisson(LES) 

■ Fractional step approach 

■ Third-order, upwind 
biased scheme for the 
convection terms 

■ Second order central 
scheme for the diffusion 
terms 

■ Second-order accurate 
pressure Poisson 
Equation 

* Samgorinsky’s scheme 
for subgrid scale stresses 


RANS CODE & TURBULENCE MODELS 

■ Control-Volume formulation to solve 
time-averaged equations 

■ Second-order accurate schemes; 
QUICK for convection, CDS for diffusion 

■ Pressure and pressure-correction 
schemes as in SIMPLER 

■ Line-by-line Gauss -Seidel Relaxation 

■ Standard k-e model; Non-linear k- e 
model; Functionalized k- e model 
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INS-3D-DNS 

■ Governing Equations: 

da dtf-j,) dCf-l) d(g-l) . 

dt dl] d£ 

dp^ q [~ 7 £) d(Y_ / J) d(W i J) 

dr dr] + d( 


/_ + 


Discretised Equation has the form: 

nAY +U l j 

At 


A 


1.5D n+l m -2/H 


+0.5D 


n-l 


■ Used fifth order upwind-biased discretization for convection; 
fourth order CDS for diffusion; Gauss-Seidal Line Relaxation 
SCheme or GMRES solver 


INS-Poisson-LES 

■ Fractional step Approach 


p"* l V - p"V” 

At 

V. — = v.v* 
VP J 

p" +1 y" +I - p s F* _ 
Ar 


C(V) + Z>(V) 


-Vp 


■Temporal accuracy increased by using either 
an Adams-Bashforth scheme, or by a second 
loop-through the steps. 
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LES: Subgrid Models 


Filtered equations of motion 


dU ; 


dt dx s 


^ (¥;) : 


\dp d 


pdx. dx 


*"T-(*,v) + V 


dxdx 


(«) 


■ Subgrid scale stresses: x » -w-m 

■ Smagorinsky’s model: Eddy viscosity 

■ Dynamic subgrid scale model: test filter 


RANS: Turbulence Models 

■ Standard k-e model 

■ Non-linear k- e model of Speziale 

■ Functionalized model of Dutta and 
Acharya 

■ Other models: RNG, RS, Durbin 
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Preliminary Flow Results 

■ Ajersch et al’s test case 



Uapp=2.75 m/s; Vjet=5.5 m/s; d=1 2.7mm; L=4d 

■ Results obtained on a 140x40x40 non- 
uniform mesh; 80 points in the hole region; 
Resolves large scale features; No subgrid 
stress models 


PARTICLE TRACES 
UjetAJo** V* 
jet centerline 


0.000 MICH 
0.OOOCG HLWtf 
4.70x10«*3 Re 
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PARTICLE TRACES 


Ujet/Uo s 1/2 
Jat centerllna 


0.000 

MRCH 

0.00 DEG 

BLPHF 
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Ra 

5.5 

TIKE 

198G5i48 

GRID 



PARTICLE TRACES 
Ujat/Uo * 1/2 
Jat cantarlina 


8.880 HACK 

0.808EC 8LPRF 

Ka 

5.5 TIKE 
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ONGOING WORK 

■ DNS on a finer grid for Ajersch’s data 

■ LES with Smagorinsky’s model 

■ RANS with improved models 
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HEAT TRANSFER IN FILM-COOLED TURBINE BLADES 


Vijay Garg 
AYT Corporation 
NASA Lewis Research Center 
Cleveland, Ohio 


MOTIVATION 


• TENDENCY TO USE HIGHER TURBINE INLET TEMPERATURE 
- TO INCREASE SPECIFIC POWER OUTPUT 
& THERMAL EFFICIENCY 


• NEED THEREFORE TO COOL THE BLADES AND SIDE WALLS 
-TO ENSURE REASONABLE LIFETIME OF MACHINE 


OBJECTIVE 


• TO DEVELOP TOOLS & MODELS TO ACCURATELY PREDICT 
HOT GAS SIDE HEAT TRANSFER TO A FILM-COOLED 
TURBINE BLADE 


• TO MINIMIZE THE INJECTED MASS FLOW FOR A GIVEN 
COOLING EFFICIENCY 


• TO DECREASE AERODYNAMIC LOSSES INDUCED BY THE 
JET-SECONDARY FLOW INTERACTION 
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EQBM.ULAT1QN 


• 3-D COMPRESSIBLE NAVIER-STOKES EQUATIONS 

• BALDWIN-LOMAX TURBULENCE MODEL WITH MAYLE’S 
TRANSITION CRITERION; ALSO TWO-EQUATION MODELS 
(k-e, q-cD, k-to); (y+ < 1 near wall) 

• BODY-FITTED COORDINATE SYSTEM; C-GRID USED 

• APPROPRIATE BOUNDARY CONDITIONS SPECIFIED 

• FOREST’S MODEL FOR LEADING EDGE HEAT TRANSFER 
AUGMENTATION DUE TO FREE-STREAM TURBULENCE 

• CRAWFORD’S MODEL FOR EDDY VISCOSITY AUGMEN- 
TATION DUE TO FILM COOLING 


FI L M COO I IRG-C-Q-NSIP ERA T I Q NS 

• EXIT AREA OF EACH HOLE DISTRIBUTED OVER SEVERAL 
CONTROL VOLUMES 


• ABILITY TO ANALYZE DIFFERENT HOLE SHAPES 


• DIFFERENT VELOCITY AND TEMPERATURE PROFILES 
SPECIFIED AT THE HOLE EXIT 


• SPECIFIED HEAT FLUX OR VARIABLE TEMPERATURE 
ON THE BLADE SURFACE 
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NUMERICAL SOLUTION 

• TRAF3D CODE (by Andrea Arnone) MODIFIED 

• CELL-CENTERED FINITE VOLUME APPROACH 

• FOUR-STAGE RUNGE-KUTTA TIME STEPPING SCHEME 

• 3-D EIGENVALUE SCALING OF ARTIFICIAL DISSIPATION 

• VARIABLE COEFFICIENT IMPLICIT RESIDUAL SMOOTHING 

• FULL APPROXIMATION STORAGE MULTI-GRID METHOD 

• EFFECT OF COOLANT VELOCITY AND TEMPERATURE 
DISTRIBUTION AT THE HOLE EXIT ON THE HEAT TRANSFER 
COEFFICIENT ON VARIOUS TURBINE BLADES 

• COMPARISON WITH EXPERIMENTAL DATA FOR ALL BLADES 
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ri or C 

Different velocity and temperature profiles at the hole exit 


2 ROWS OF HOLES 
35° ANGLE CHORD WISE 
90° ANGLE SPANWISE 
3d SPANWISE SPACING 



5 ROWS OF STAGGERED HOLES 
90° ANGLE CHORDWISE 
45° ANGLE SPANWISE 
7J>d SPANWISE SPACING 


HOLE DIA. d = 0.99 mm 
4d ROW SPACING 


C3X vane and cooling hole details 
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I I I I I I I I 


1 I I I I I I I I 


h = 1 1 35.6 W/m 2 -K a experimental data 


15 [_ 281x45x81 grid 


polynomial profile 
1/7th power law 


k & /\/ h 




-0.5 

Pressure 


0.5 

Suction 


Effect of Velocity & temperature profiles at the hole exit 
(C3X vane: B-L turbulence model; case 44355) 
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VKI rotor and cooling hole details 
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Pressure Suction 

s/c 

Effect of Velocity & temperature profiles at the hole exit 
(VKI rotor: B-L turbulence model; 2.07% coolant) 



s/c 

Effect of Velocity & temperature profiles at the hole exit 
(VKI rotor: B-L turbulence model; 3.09% coolant) 
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s 


Effect of Velocity & temperature profiles at the hole exit 
(ACE rotor: B-L turbulence model; case 6115) 


• EFFECT OF BLADE ROTATION 

• EFFECT OF DIRECTION OF COOLANT INJECTION FROM THE 
SHOWER-HEAD HOLES 
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Demonstration of Advanced Film Cooling Model to a Turbine Blade Thermal Design Problem 
(Adjusting Coolant Injection Angles to Maximize Cooling Effectiveness) 


ADIABATIC EFFECTIVENESS 

1.0 


a) 

0.5 

b) 

c) 

0.0 



Q. = 5200 rpm 
Tc/To = 0.5 
mc/mo = 3.65% 
ic/lo = 3.55% 


UTRC 


All Shower-Head Holes 
inject towards the Hub 

All Shower-Head Holes 
inject towards the Tip 

All Shower-Head Holes 
inject towards Tip except 
8 near hub on suction side 
that inject towards Hub 



Tip 


• EFFECT OF TURBULENCE MODELING: 

Baldwin-Lomax model 
Coakley’s q-co model 
Chien’s k-s model 
Wilcox’s k-co model 


• COMPARISON WITH EXPERIMENTAL DATA FOR ALL BLADES 



Pressure Suction 

S 

Various Turbulence Models: C3X vane near mid-span (case 44135) 
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Various Turbulence Models: C3X vane near mid-span (case 44155) 



Pressure Suction 

S 

Various Turbulence Models: C3X vane near mid-span (case 44355) 
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Grid on the blade surface, and C-grid on 
the hub. Suction surface is in the front. 
Alternate grid lines shown. 



Shape of film cooling holes at exit on 
the ACE rotor surface 
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Leading Edge 



Heat flux gauges and cooling hole locations on 
the projected blade surface with each of the three 
chordwise rows of gauges on a separate blade 
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Pressure Suction 












Pressure Suction 

S 



Pressure Suction 

S 



Pressure Suction 

S 

Nusselt Number on the Blade Surface for Cooled Case 72 
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Pressure Suction 


S 



Pressure Suction 


S 



Pressure Suction 

S 

Nusselt Number on the Blade Surface for Cooled Case 73 
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ALSO STUDIED 


• EFFECT OF SPANWISE PITCH OF SHOWER-HEAD HOLES 

• EFFECT OF COOLANT TO MAINSTREAM MASS FLOW RATIO 

• EFFECT OF COOLANT TO MAINSTREAM TEMPERATURE 
RATIO 


CONCL U SI ONS 

• A LARGE CHANGE IN HEAT TRANSFER COEFFICIENT AT 
THE BLADE SURFACE CAUSED BY DIFFERENT VELOCITY 
AND TEMPERATURE PROFILES AT THE HOLE EXIT 

• DIFFERENT EFFECTS OBSERVED ON THE PRESSURE AND 
SUCTION SURFACE DEPENDING UPON THE BLADE & HOLE 
SHAPE (CONICAL OR CYLINDRICAL) AND HOLE LOCATION 

• BLADE ROTATION AND DIRECTION OF INJECTION HAS A 
SIGNIFICANT EFFECT ON THE FILM COOLING EFFECTIVENESS 

• BALDWIN-LOMAX TURBULENCE MODEL SEEMS TO BE AS 
GOOD AS ANY TWO-EQUATION MODEL, AND IS 40-50% LESS 
EXPENSIVE 
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FUTURE DIRECTIONS 


• DETAILS OF IN-HOLE AND NEAR-HOLE PHYSICS REQUIRED 
UNDER CONDITIONS RELEVANT TO THE GAS TURBINE 
(DENSITY RATIO; CURVATURE EFFECTS; HIGH Tu; etc.) 

• EXPERIMENTAL DATA ON ROTATING FILM-COOLED BLADES 
REQUIRED 

• INCORPORATE HOLE-EXIT PROFILES BASED ON RECENT 
EXPERIMENTAL DATA 

• USE A MULTI-BLOCK CODE TO DEVELOP HOLE-EXIT 
PROFILES ON A ROTATING BLADE 

• INCORPORATE SIMPLIFIED FILM-COOLING MODEL TO 
SPEED-UP THE CODE 
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THE EFFECT OF TABS ON A JET IN A CROSS-FLOW 


Khairul Zaman and Judith Foss Van Zante 
NASA Lewis Research Center 
Cleveland, Ohio 


Objective: 

A tab placed on the leeward side of the 
nozzle was expected to increase jet 
penetration into the cross-flow. An 
experiment at UTRC showed 
insignificant effect. The primary 
objective of the present study was to 
confirm and explain the ineffectiveness. 


Approach: 

• Experiments in a low speed wind tunnel. 

• Conduct hot-wire measurements for 
mean velocity and streamwise vorticity 
fields. 
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Results: 

• Ineffectiveness of tab, when placed on 
leeward side, is confirmed. 

• Static pressure distribution around the 
nozzle provides an explanation for the 
ineffectiveness. 

• Tab placed on the windward side of the 
nozzle is found to reduce the jet 
penetration. 

Milestones: 

• 5/97 Experiments to optimize tab 
geomery for maximum reduction in jet 
penetration. 

• 5/98 Active control for increasing jet 
penetration and mixing. 


SCHEMATIC OF ’DELTA-TAB’ 
FITTED TO CIRCULAR NOZZLE 
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Fric & Roshko JFM 1994 


Counter-rotating 
Vortex Pair 




A Tab on the downstream edge 
may be expected to reinforce 
the Counter-rotating Vortex Fair 
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Liscinsky, True & Holdeman AIM 95-2998 
(Experiment at UTRC) 
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(Pst - p A )/VipU T J 


STATIC PRESSURE DISTRIBUTION NEAR THE NOZZLE EXIT 
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p dz dy 
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POSSIBLE INVESTIGATION FOR FILM COOLING 

APPLICATION 


THE FACT THAT A TAB PLACED ON THE WINDWARD 
SIDE REDUCES JET PENETRATION IS OF INTEREST IN FILM 
COOLING APPLICATION. WE WOULD LIKE TO PURSUE THIS 
FURTHER THROUGH EXPERIMENTS INCLUDING: 

(1) BASIC STUDIES ON TAB GEOMETRY OPTIMIZATION (LOW 
‘J% ANGLED JET, ARRAY OF JETS). 

(2) HEAT TRANSFER. 

(3) TESTS FOR TURBINE BLADES AND OTHER APPLICATIONS. 


190 



FILM-COOLING HEAT-TRANSFER MEASUREMENTS USING LIQUID CRYSTALS 

Steven Hippensteele 
NASA Lewis Research Center 
Cleveland, Ohio 


The Transient Liquid-Crystal Heat-Transfer Technique 

2-D Film-Cooling Heat-Transfer on an AlliedSignal Vane 

Effects of Tab Vortex Generators on Surface Heat Transfer 
Downstream of a Jet in Crossflow 


Transient Heat Transfer Technique 


Based on 1-D heat conduction out of a semi-infinite wall: 
0 = \-e^ erf eft 

T -T 

d = — — - J£ - 
T -T 



h4t 

■Jpck 


erfc fi = complimentary error function 

i=initial, LC=liquid crystal, R=recovery 

p=wall density, c=specific heat, k=thermal conductivity 
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Transient Liquid-Crystal Heat-Transfer Technique 


Vane surface coated with thermochromic liquid crystals 

Initially heated in oven to a uniform temperature and cooled during the test (see oven 
sketch and photo) 

Isotherms (revealed by the liquid crystals) are video taped and used to solve the 
1-D conduction problem in the semi-infinite test surface 

By operating the coolant and air at a constant room temperature, only the 
hydrodynamic effects will be included the heat transfer coefficients 

By varying the coolant temperature, the cooling effectiveness can be determined 
(ref. Vedula & Metzger) 

Results: high-resolution heat-transfer and film-cooling effectiveness maps of the 
surface downstream of the film cooling holes 


Test Section Oven 
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2-B Film-Cooling Heat-Transfer on am 
AlliedSignal Turbine Vane 


The unwrapped midspan geometry of the AlliedSignal TFE I042-70A high pressure turbine 
first stage vane has been offered for detailed heat transfer studies and Lewis, This vane 
contains 12 rows of film cooling holes. A linear cascade will be fabricated using this vane 
geometry for the experimental study of surface heat transfer using the transient liquid-crystal 
heat-transfer technique. The vane size will be about 4 times actual size so that both the exit 
Mach number and the exit Reynolds number can be maintained. The coolant will be SF 6 so 
that the coolant-to-gas density ratio can be maintained and changed. The test vane surface will 
be coated with thermochromic liquid crystals and initially heated in an oven. The isotherms 
revealed by the liquid crystals, that result from the cooling during the test period, are video 
taped and used to solve the 1-D conduction problem in the semi-infinite test surface. By 
operating the coolant and air at a constant room temperature, only the hydrodynamic effects 
will be included the heat transfer coefficients. By varying the coolant temperature, the cooling 
effectiveness can be determined. From this comes high-resolution heat transfer and film 
cooling effectiveness maps of the surface downstream of the film cooling holes. 
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2-D Film-Cooling Heat-Transfer on an 

AlliedSignal Vane 


Detailed heat transfer studies to validate film cooling HT 
models (Garg) 

Using the transient liquid-crystal heat-transfer technique 


The Vane Cascade 

AlliedSignal offered the unwrapped midspan geometry of their TFE 
1042-70A high pressure turbine first-stage vane for study and non- 
proprietary reporting (Winstanley) 

Contains 12 rows of film cooling holes 

4 times actual size stereo-lithographic model (see photo) 

Real exit Mach number and Reynolds number will be tested 

Coolant will be SF 6 so real coolant-to-gas density ratio can be tested 

A linear cascade of vanes will be used for the experimental study (heat 
transfer, cooling effectiveness, hot-wire surveys, & velocity surveys with 
high free-stream turbulence (Ames)) 
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considerable effect on the downstream flow (Zaman and Foss). The tab-generated vorticity 
reduces the penetration of the jet which may be significant in film cooling applications. For 
this reason experimental measurements have been initiated using the transient liquid-crystal 
heat-transfer technique in which an acrylic test surface containing an air jet coated with 

r 
■& 



Effects of Tab Vortex Generators on Heat 
Transfer Downstream of a Jet in Crossflow 

Tab vortex generators (windward side of a jet in a crossflow) have shown 
considerable aerodynamic effects on the downstream flow (Zaman and Foss) 

The tab-generated vorticity reduces jet penetration which may be significant in 
film cooling applications 

Therefore, experimental measurements have started using the transient liquid- 
crystal heat-transfer technique (see facility sketch and photo) 

High-resolution heat-transfer maps will be made for the surface downstream of 
the jet (and eventually for film cooling holes) (see 2 isothermal/h maps) 

Initial data, as available (shown in video) 


Transient Liquid-Crystal Heat Transfer Facility 
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Jet With Tab 


Preliminary Data 




TURBINE VANE HEAT TRANSFER, FILM COOLING, AND TURBULENCE EXPERIMENTS 

Forrest E. Ames 
Allison Engine Company 
Indianapolis, Indiana 


PRESENTATION OVERVIEW 


I. Introduction 

o Motivation 
o Objectives 

II. Experimental Apparatus and Approach 

III. Experimental Results 

o Suction Surface Film Cooling Comparisons 
o Pressure Surface Film Cooling and Velocity Profiles 
o Suction Surface Heat Transfer with Film Cooling 
o Pressure Surface Heat Transfer with Film Cooling 

IV. Conclusions and Thoughts for Future Work 
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INTRODUCTION 


I. Several years ago a test was run at Allison on an advanced film 
cooled vane design. 

o The tests indicated we were overpredicting film cooling 
protection on the pressure side, 
o Results suggested we were underpredicting surface heat 
transfer downstream from the pressure surface film cooling. 

II. The current study was undertaken to improve our understanding of 
the film cooled pressure surfaces with high inlet turbulence. 

III. Study’s objectives were: 

o Examine influence of high inlet turbulence on vane film cooling 
o Document local heat transfer rate and BL with jet injection 
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EXPERIMENTAL APPROACH 



Schematic of four vane C3X cascade 
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.635 


14.60 



II 

Cross Sectional View 


Side Panel 


Back Pane! 


Schematic of mock combustor turbulence generator 


Suction surface film cooling holes 



cooling array 

C3X film cooled vane thermocouple positions 
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FILM COOLING RESULTS 

o Comparison between 
L’Ecuyer correlation and 
data show good agreement 
both in level and in trends. 

o Agreement both suggests 
suction surface film cool- 
ing is comparable to flat 
plate data and gives 
confidence in experimental 
method. 

o Data show higher near hole 
effectiveness levels for low 
velocity ratios but lower 
protection downstream. 
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■Had, (Taw-Tr)/(Tco-Tr) 



o Typical range of velocity ratios for first stage vane suction surface is quite 
narrow due to low Pc/Pt, low local Ps, high Tg/Tco, and hole discharge loss. 
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o Very high acceleration 
appears to cause large rise 
in near hole r| (VR= 5, Tu= 01) 
and significant reduction in 
the far field. 

o Data show large reduction 
in rj with Tu 

o Near hole data have a large 
uncertainty due to corrected 
rj based on 2-D conduction 
analysis. Low Tu condition 
for low VR has largest cr. 

o Conduction correction for 
X/D > 8 is very small. 
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■Had, (Taw-Tr)/(Tco-Tr) 



0 10 20 30 40 50 60 70 80 

X/D 

Pressure surface film cooling, low & high turbulence, single row, VR = 0.5 

o Overall result of combination of high turbulence level and high acceleration 
is very low effectiveness level on the PS. 
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0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Y(cm) 

Low turbulence PS velocity profile, VR = 0.5, one row 

o Low turbulence PS boundary layer shows a large 
velocity deficit downstream from jet but little effect 
on the profiles taken at Z/D’s of 0.75 and 1 .5. 

o Turbulence level in jet moderate in intensity (6%) 
peaks in shear layer (7%) and low in between. 

o Dissipation moderate in jet very high in shear layer. 
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— pi 

cbv5a2, VR=0.5, Z/D=0„ Tu=.12, Lu=3.36cm 
Cbv5q2, VR=0.5, Z/D=,75, Tu=.12, Lu=3.36cm 
cbv5m2, VR=0.5, Z/D=1 .5. Tu=.12, Lu=3.36cm 
cbvOal, Tu=.12, Lu=3.36, BaseVane 



m 


0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Y (cm) 


High turbulence PS velocity profile, VR = 0.5, one row 

o Time averaged high turbulence PS boundary layer 
shows excellent uniformity across entire span by 
anX/DoflO. 

o Turbulence level is uniformly high for all spanwise 
conditions 

o Dissipation level 50% higher than basevane at .1 cm 
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o Data for two rows of holes 
on the PS at low and high 
turbulence are compared to 
One row data superposed at 
a velocity ratio of 0.5 . 

o Comparison of low and high 
turbulence data show a 
substantial reduction in rj 
due to turbulence. 

o Double row and superposed 
single row comparison for 
high Tu indicates no benefit 
from staggered row due to 
extra blockage effect. 
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o Similar to the single row data, combination of turbulence mixing and high 
acceleration produces very low effectiveness level on PS compared to SS. 



o Data for two rows of holes 
on the PS at higher VRs 
show higher rj values than 
superposed values from 
one row. 

o For low turbulence, higher 
blockage from the double 
staggered row substantially 
reduces penetration of the 
downstream jet for VR = 1 .5 

o Velocity profiles also show 
less liftoff for the upstream 
jet. 
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T|ad, (Taw-Tr)/(Tco-Tr) 



0 10 20 30 m 40 50 60 70 80 

Pressure surface effectiveness, 2 rows, 30° holes, VR = 1.5, DR=0.94 

o While the staggered row arrangement produces a substantial enhancement 
to film cooling on the PS, turbulent mixing from high local turbulence levels 
produce substantial reductions in effectiveness. 
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o The data for two rows can 
be compared for low and 
high turbulence by ratioing 
r| for similar conditions. 

o The average ratio of the four 
conditions decreases with 
X/D indicating how turbulent 
mixing gradually but 
increasingly reduces x\. 

o Generally, all the ratios of r\ 
decrease with increasing 
X/D. Near hole ratios for a 
VR of one show higher r\ 
values at an X/D of four for 
the high turbulence case. 
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Ratio of high and low turbulence effectiveness, PS, 2 rows, 30° holes, DR=0.94 

o High turbulence t| values may benefit from enhanced turbulent mixing in the 
near hole region. values in the near hole region are also more uncertain 
due to conduction effects. 
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HEAT TRANSFER RESULTS 


o Two staggered rows of 
film cooling produce an 
average augmentation peak 
of about 18 percent above 
the turbulent baseline. 

o The heat transfer quickly 
decays to the base level. 

o The low turbulence case 
produced a comparable 
level of augmentation. 

o A single row of holes 
produced an augmentation 
ratio of about 1 .1 1 for both 
high and low turbulence 
levels. 
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Stanton number 



-15 -10 -5 0 5 10 15 20 

Surface distance (cm) 

Suction surface heat transfer, 2 rows, 30° holes, high turbulence 

o The level of augmentation was found to be reasonably consistent over the 
range of velocity ratios investigated for either one or two rows. The heat 
transfer increase for the showerhead array was also low past transition. 
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o In spite of an augmentation 
level of 60 percent above 
the low turbulence case, the 
heat transfer was increased 
by up to an additional 46 
percent over the base case 
for a VR of 1.5. 

o Generally, heat transfer 
augmentation ratio scaled 
on velocity ratio. 

o Similar to the augmentation 
on the suction surface 
boundary layer, the effect of 
the jet quickly dissipated in 
the downstream direction. 
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Stanton number 
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Pressure surface heat transfer, 1 row, 30° holes, high turbulence 


o The peak augmentation on the pressure surface for the single row at a 
velocity ratio of 0.5 was about 16 percent compared to a peak of 1 1 percent 
on the suction surface for a single row and a comparable velocity ratio. 
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o The augmentation curves 
for velocity ratios of 1 .0 and 
1 .5 were very similar in 
form. 

o The augmentation curve for 
a velocity ratio of 0.5 had a 
delayed peak and a slower 
decay. 

o The delayed peak level and 
decay were believed due to 
the shear layer that formed 
and caused an increased 
production of turbulence. 

o At an X/D of 9, dissipation 
rates for velocity ratios of 
0.5 and 1.0 were close. 
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st/sto 


1.5 



Pressure surface heat transfer, 1 row, 30° holes, high turbulence 

o The dissipation level at Y of 0.1 cm for a VR of 1 .5 was about three times the 
level determined for a VR of 1 .0 and about 3.5 times higher than for a VR of 
0.5 indicating no simple relationship between dissipation and augmentation. 
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HEAT TRANSFER RESULTS 


o Two rows of holes produce 
about twice the increase 
in peak heat transfer as a 
single row. 

o The augmentation curves 
look quite similar with the 
curves for a single row. 

o The decay rate for a double 
row is quicker than the rate 
of decay in augmentation 
for a single row. 

o The augmentation curve for 
the velocity ratio of 0.5 has 
a different shape and decay 
rate than the higher VR data 
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st/sto 



X/D 


Pressure surface heat transfer, 2 rows, 30° holes, high turbulence 

o The dissipation levels at Y of 0.1 cm for VR’s of 0.5 and 1.0 are noticeably 
higher than single row values as measured at X/D of 9. The dissipation 
levels determined for the VR of 1 .5 are comparable to single row data. 
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CONCLUSIONS 


o For this C3X vane, suction surface film 
cooling produced levels of rj similar to a 
flat plate configuration for the low Tu case. 

o High pressure surface acceleration raised 
near hole rj over suction surface but r\ level 
diminished with continued velocity increase. 

o High relative level of turbulence on PS 
caused rapid and continuing depletion of 
film cooling coverage. 

o Two staggered rows of film cooling produced 
substantial enhancements over a single row 
for high VR’s. 

o Peak heat transfer augmentation for suction 
surface was about 1 1 percent for one row and 
18 percent for two for both turbulence levels. 

o PS augmentation scaled on VR and number of 
rows with peak augmentation levels of up to 95 
percent over high turbulence base for two 
staggered rows at a VR of 1.5. 
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o Flows with strong continuous acceleration 
need to be studied at various Tu levels, 
(different positions along a vane PS) 

o Effect of relative hole size on augmentation 
also needs to be evaluated. 
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INVESTIGATION OF ROTOR WAKE EFFECTS ON FILM COOLING 

James D. Heidmann 
NASA Lewis Research Center 
Cleveland, Ohio 

Ou t li n e 

- Motivation and goals 

- Experimental study 

- Computational study 

- Wake-affected film cooling model 

- Conclusions and recommendations 

Motivation and Goals 

- Effect of wake passing on film cooling is not well understood 

- Lack of understanding may result in excess coolant usage 

- Focus on showerhead cooling due to sensitivity to wake passing 

- Temporal, spanwise, and chordwise resolution of test data 

- Seeking physics-based model to account for wake passing effect 

- Model should be applicable over a broad range of parameters 
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Best-Fit Equations for Air and CP 2 


Air: 

T) = 


1.0 + ( 


0.761 

0.054x 


MS(1.0 ± 0.139M ± 0.286St) 


) 


0.792 ± 0.033 


Ati^ = 0.0068 


C0 2 : 
T1 = - 


0.948 


1.0 + ( 


0.094x 


MS(1.0 ± 0.241M ± 0.144St) 


) 


0.762 ± 0.014 


Arinas = 0.0074 


0.093St 


0.095St 
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Film effectiveness Film effectiveness 


Correlated Snan-Avera 




Distance from leading edge, x/d 
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Film effectiveness Film effectiveness 




Distance from leading edge, x/d 







Conclusions and Recommendations 


- Primary wake effect is linear reduction in T| with St 

- Secondary wake effect is skewing of suction/pressure side cooling 

- Steady computations match experimental Nu, but overpredict r) 

- Unsteady computations elucidate wake/film interaction 

- Model may be used to estimate wake passing effect 

- Need boundary layer and full stage experiments 

- Need resolved film hole and full stage unsteady computations 

- Need validated turbulence models for film cooling 
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UNSTEADY HIGH TURBULENCE EFFECT ON TURBINE BLADE FILM COOLING HEAT 
TRANSFER PERFORMANCE USING A TRANSIENT LIQUID CRYSTAL TECHNIQUE 


H. Du, S.V. Ekkad, and J.C. Han 
Texas A&M University 
College Station, Texas 


OBJECTIVES 


To study detailed distributions of heat transfer coefficients and film effectiveness 
for a gas turbine blade, including film cooling hole region, under the effects of 
upstream flow conditions: 

• Free stream turbulence 

• Unsteady wakes 

• Unsteady wakes with trailing edge coolant ejections 

and the effects of film cooling injections with air and C0 2 as coolants: 

• Blowing ratios 

• Coolant-to-mainstream density ratios 

• Advanced film hole shapes. 
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Sketch of Test Section 



Camera Arrangement 
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Axial 


Liquid Crystal 
Coated Blade 
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APPROACH 


• A transient liquid crystal technique has been developed for this study. 

• In the experiments, the liquid crystal coated blade is heated to a uniform 
temperature, and then suddenly exposed to a pre-set mainstream flow condition. 

• A Color Frame Grabber system, through the RGB cameras, is used to record 
the time of the color change of the liquid crystal on the blade surface. 

• The color change times are related to heat transfer coefficients and film 
effectiveness based on the solution of a 1-D semi-infinite transient conduction 
model with a convective boundary conditions. 


THEORY 


ONE-DIMENSIONAL SEMI-INFINITE SOLID ASSUMPTION: 





Initial condition: 

at t=0, T=T l 


Boundary conditions: 

at x=0, -k— = h(T w -TJ; as T=T { 

dx 

Solution at x=0: 


T w~ T i , ,h 2 at . . .hsfca 

* 1 -exp( — —)erfc(— — ) 

T -T, k 2 k 

m l a 
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FILM COOLING CONCEPT 


Mainstream 

u m’Tm 



Coolant U c , T c 


q // 0 - h g (T m -TJ with no film cooling 


q y/ - h(T f -TJ with film cooling 


T -T 

7 7 = — — L, film effectiveness 
T -T 

m c 


r!L - h T r T ~ 

q " 0 h o T m' T w 
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For Film Cooling Case: T m = T r 
Solution: 


T w T i , ,h 2 at , r , 

- 1 -exp( )erfc( 

T f~ T i k 2 


where T f = local film temperature 


it 


T f * V(T C ~TJ+T m = TfT e m-rj)T m 
t|— local film effectiveness 

Solution: 






r^-r, = /2 -exp(^)er/c(^L)l[rjT c+ (l -Tj)T m -TJ 


Measured Region on Liquid Crystal Coated Blade 


Thin Foil liquid Crystal 

Heater Coated Regoin 
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Heater Box Cross Section 


Heater Box 




Conceptual View of Unsteady Wake Effect on Blade 
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\7 (m/s) 



0 5 10 15 20 0 5 10 15 20 


Time (ms) Time (ms) 

Local-to-exit velocity ratio (V/V 2 ) distributions on the test blade, 
profiles for V(t), V . Tu , under wake effect at cascade inlet 
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Effect of Unsteady Wake on Spanwise- Averaged Nusselt Number 
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Effect of Unsteady Wake on Spanwise- Averaged Film Effectiveness 
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Blowing Ratio Effect on Spanwise- Averaged Nusselt Number at Given 
Wake Condition 
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Density Ratio Effect on Spanwise-Averaged Film Effectiveness at Given 
Wake Condition 
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Conceptual View of Effect of Unsteady Wake with Trailing Edge Ejection 



Sketch of Rotating Rod with Ejection Holes 


Rotating Direction 




Rotor 


256 




V(m/») 









X/PL X/SL 


Combined Effect of Unsteady Wake, Free-Stream Turbulence and 
Coolant Ejection on Spanwise- Averaged Nusselt Number 
for Re = 7.6 x 10 5 


MILESTONE 


The effects of upstream flow conditions and film cooling injections have been 
studied. The tasks for the coming year are: 

• To design and fabricate a test blade with several rows of shaped film holes. 
These advanced film hole shape designs will be decided after consulting with 
GE for the optimum shapes. 

• To measure detailed heat transfer coefficients and film effectiveness for the 
blade with shaped film cooling holes under the combined effects of free-stream 
turbulence, unsteady wake and unsteady wake with coolant ejection. 
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RESULTS 


• Nusselt numbers increase with an increase in free-stream turbulence level, and 
also increase with the addition of unsteady wakes. Adding grid-generated 
turbulence to unsteady wakes further increases Nusselt numbers on the 
downstream blade surface. 

• The trailing edge ejection is to increase Nusselt numbers on the front parts of 
both the blade suction and pressure surfaces. 

• Film injection promotes earlier laminar to turbulent boundary layer transition 
on the suction surface and also greatly enhances local Nusselt numbers. 

• Unsteady wakes slightly enhance Nusselt numbers but significantly reduce film 
effectiveness on a film-cooled blade compared with a film-cooled blade without 
unsteady wakes. 

• An increase in blowing ratio increases Nusselt numbers for both coolants. Film 
effectiveness is highest at a blowing ratio of 0.8 for air injection and at a 
blowing ratio of 1.2 for C0 2 injection. 
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GROUP DISCUSSION - December 12, 1996 

SESSION II. FILM COOLING 


Rav Gaueler - Facilitator Nirm Nirmalan - Scribe 

We have new detail computational methods. That is not our end product. Vision in my 
mind is to have something that can be used by designers. Final Goal: we are trying to 
find out what is going on. NASA Role & Vision. 

Terry Simon 

Are shaped holes all product specific? Are they better left inside the company? 

Rav Gaueler 

NASA, under AST, is working with an AlliedSignal vane that has shaped holes. 

Terry Simon 

Who decided the shape? 

Rav Gaueler 
AlliedSignal did. 

Yong Klim 

From A/F contract, the shaped hole. Film cooling (FC) hole shape depends where on the vane it 
is to be located. 

Rav Gaugler 

Lots of shaped holes in upcoming papers in the IGTI. It is product specific in the AST program. 
Bob Bergholz 

AST has both product and non-product specific. Testing will occur and will be shared within the 
AST community. 
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Phil Ligrani 

Free stream turbulence & wake passing } Work is being done. 


Shock Interactions 
Flow Interactions 
Relative Importance 

Shock Wave 

Shaped Holes 

Luzene Zhang 


Anything going? 

Oxford: Martin Oldfield 

We are interested, but it is difficult. 

Need to study with curvature. 


Shock waves. Interested but they are very difficult to study. 

Also in shaped holes, you need to consider acceleration and curvature influences. 

Rav Gaugler 

There are shaped hole data. 


J. C. Han 

We are trying to find optimum shaped holes. P&W has the most. 

Boris Glezer 

What are the shape of holes? One good for one airfoil is not the same for another airfoil. Need 
to address this effect of curvature. 

Jim Heidmann 

Is it the curvature effect along with the pressure gradient effect? 

Ron Bunker 

On the AlliedSignal blade the shaped hole looks almost standard. Pretty good! 

Younes Makki 

P&W (maybe) showed improved effectiveness. 

Joel Wagner 

Testing in Florida. All kinds of shaped holes. Manufacturing (EDM) is very expensive. 
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Younes Makki 


Recast and EDM are expensive. 
Bob Simoneau 


Recall of Steve Papell cusp shaped FC holes to produce vorticity. P&W took it and got data. 
They looked at it and other shaped holes. Looked good. 

Mine-King Chvu 

Manufacturing is an issue. Variability of shaped holes are vast. Leylak has done some 
modeling. Need to stay with a hole. Should standardize their shape for study. 

Luzene Zhang 

Lot of data in the 70’s. Shaped holes 1) reduces the exit blowing ratio (BR), 2) improves 
spanwise spreadability, and 3) improves exit vorticity and therefore the stickability. Our test 
shows no difference between shaped and round holes. At high BR there is a difference. 

Bob Simoneau 


[To Terry Simon] Are you planning to look at FC inlet air at 90 degrees to the hole? 
Terry Simon 


Yes, we have that test section installed now. The flow startes to turn streamwise even before it 
leaves the FC hole. Taking data in the hole is hard and using a hot wire is uncertain. May need 
to use LDA. We would like to work on compound and shaped holes. 

Yong Kim 

[To Terry Simon] Flow was measured but pressure drop (CD) was not. 

Terry Simon 


We did not measure pressure drop. We metered the flow. Velocity is not too high. 

[To Karen Thole] You were involved in an experiment where you had a recirculating flow? 

Karen Thole 


CDs were different. Tests were done in Germany. Hole shape does not help diffusing very 
much. Inlet direction does matter in shaped holes (CFD). 
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Terry Simon 
We should do CD. 

Yong Kim 

Shaped holes are valuable. 

Karen Thole 

Measured static pressure around the hole: there is variation. 

Terry Simon 
Thanks. 

J. C. Han 

Can anyone comment on how do you compare shaped and compound holes? 

Luzeng Zhang 

Bogard has a paper on comparing holes. 

Karen Thole 

Bogard did compound hole. 

Khairul Zaman 

Measuring pressure drop in the tunnel is too difficult because it is too small. Need to use an 
orifice to meter the flow. 

Terrv Simon 

That is what we do. But I need to measure static pressure for CD. Measure from inlet total 
(plenum) pressure to external static pressure. 

Chi-Rong Wang 

Upstream recirculation effects on incoming flow. 

J. C. Han 

We have some data. No data to compare with. Could be a problem. 


264 



Mine-King Chvu 

Liquid crystal heat transfer is not good near film cooling holes: ID conduction does not work. 
Check on data reduction method. 

Joel Wagner 

3D conduction is going on in and around the film cooling hole. No work done yet. 

Luzeng Zhang 

Goldstein has a paper. 1993 Turbine Conference (ASME , Orlando)? 

Ming-King Chvu 

Liquid crystal is a work horse. Napthalene will change shape. 

J. C. Han 

ID conducton will take care of it. 

Simon 

We did not see recirculation zone with a 30 degree angled hole. Acharya was able to see it 
because it was a 90 degree angled hole. 

Unknown 

Literature on jets and nozzles. There are shapes. Any consideration of this. 

Khairul Zaman 

Jets spread but film cooling holes should not. 

Simon 

Forrest Ames showed mixing improved the sticking effect of film cooling flow 

Forrest Ames 

No. 

Yong Kim 

Liquid crystal heat transfer: if temp, is periodic, the solution is wrong. I had a problem with it. 
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Satish Ramadhvani 

Plenums in experiments are smooth. Will it have trips and augmentors in it? 

Rav Gaueler 
Any answers? 

Bob Bereholz 
Maybe. 

Sai Lau 

Film cooling, rotation? 

Luzeng Zhang 
Mitsubishi started work. 

Ed North 

Lots of problems with it (rotation data). Rotation data is vague and difficult, rare. 

Sai Lau 

What are general shapes? Flat plate? Rotation and buoyancy effects? 

Boris Glezer 

Avoid film cooling in rotating parts. 

Joel Wagner 

We studied with SF 6 . Blade, Mike Blair: density, buoyancy effects. On pressure side, strong 
effects. Suction side do not remember. Had a single hole. 

Robert Norton 

In a low speed rig, is it applicable to a real engine, like A/F, MIT or Mike Dunn? 

Joel Wagner 

Dust conditions give good flow visualization. 
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Satish Ramadhvani 

What about studying physics. A rotating cone: Coriolis, buoyancy? Will that help? 

Rav Gaueler 

[To Vijay Garg] Can you simulate density difference? 

Viiav Garg 

Can put in temperature to simulate density difference. ACE rotor: some holes were choked. 
When CFD group looks at data, lots of questions. For 2 equation model: Ali Ameri does not 
know what is coming out of holes. 

Rob Norton 

Why not do a parametric study to see the effects. 

Viiav Garg 

It only makes a difference very near the holes: have some effects (Tu & Scale), but there is no 
data very near the holes. 

Bob Simoneau 

Comments about Satish Ramadhyani’s cone idea: will people find it useful? Is it an academic 
exam? 

Joel Wagner 

That may be the way to get the effect of curvature. 

Boris Glezer 
Maybe sometimes. 

Rav Gaueler 
Adjourn. 
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Session HI 


Coolant Flow Integration & Optimization 




MULTIDISCIPLINARY OPTIMIZATION OF COOLED TURBINE DESIGN USING 

3-D CFD ANALYSIS* 


Aditi Chattopdhyay, John N. Rajadas, 
and Shashishekara Talya S. 

Arizona State University 
Tempe, Arizona 

OBJECTIVES 

Minimize turbine blade temperature and optimize 
coolant flow management using Multiobjective, 
Multidisciplinary Optimization Procedure 

Design parameters used in optimization 

* Blade shape 

* Coolant path geometry 

* Film cooling parameters 


DESIGN VARIABLES 

Blade chord length 
Surface geometry parameters 

- Number of splines 
Spline parameters 

Coolant path configuration 

- Coolant hole geometry parameters 

* Presented at the conference by Yong Kim, AlliedSignal Aerospace. 
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DESIGN VARIABLES 
(contd) 

Film cooling parameters 

- Position of holes on blade surface 

- Shape and size of holes 

- Angle(s) of injection 

- Injection rate 


OPTIMIZATION FORMULATION 

Minimize 

Maximum blade temperature (T max ) 

Subject to constraints on 

Aerodynamic performance 
Average temperature (T ave £ T aH ) 

Blade thickness 
Trailing edge angle 
Throat width 
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BLADE MODEL 


Interior 


Coolant Paths 



ANALYSIS 

CFDfor flow field 

- RVC3D / RVCQ3D (NASA Lewis Research Center) 

Heat transfer 

- Finite element analysis in blade interior (ASU) 

Multiobjective, Multidisciplinary Optimization 
Procedure (ASU) 

- Kreisselmeier - Steinhauser (K-S) function approach 


Multilevel format 




CFD/EXPERIMENTS TO OPTIMIZE FILM-COOLING HOLE SHAPE DESIGN TOOLS 


YongKim 

AlliedSignal Aerospace 
Phoenix, Arizona 


Advanced Subsonic Technology 

Area of Interest: 7 Affordable Turbine Cooling 
Technology 


Goal: 

• Reduce DOC by 3 to 5% by Increased TRIT of 
1 00 to 200 Degrees F Without Excessive 
Increase in Cooling Flow 

Approach: 

• Produce Benchmark 3-D Unsteady Film 
Cooling Heat Transfer Measurements 

• Validate 3-D CFD Models to Benchmark Data. 
Increase Film Cooling Effectiveness with 
Advanced Hole Shape 

• Initiate Multidisciplinary Optimization (MDO) 
Applicable to Airfoil Cooling 

Milestones: 

• Complete Benchmark Measurements 3Q99 

• Complete Film Cooling 3-D CFD 

Model Validation 4Q99 

• Complete Initial MDO of Turbine 

Airfoils 4Q99 





5 220-02 
5 : 23Z+-02 
5 1 Sc-02 
5 16S-02 
S.l4e-02 


5 . 120*02 















Year 1 

Year 2 

Year 3 





Task 1 

Development of Multidisciplinary Design 
Optimization Procedure 


Task 2 

Turbine Design Optimization (2-D) 


Task 3 

3-D Aerodynamic Optimization (inviscid) 


Task 4 I 

Extension of Task 3 to H 
viscous effects D 


CFD AS A FILM-COOLING DESIGN TOOL 



D 

1 supply 

• Location (Re Effect) \ 

• Hole Size/Shape 

• Discharge Coefficient J 

• Curvature V 

• Density Ratio 
•Tu Effect 

• Surface Roughness 

• Unsteadiness 

• Heat Transfer Coeff. J 


Pressure Ratio 

T 

Mass Flux Ratio 

5 

Evaluate Heat Flux 

I 


Thermal Analysis 


i 

r 

Adjust 
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APPROACH 



CODE RANKING/CAUBRATION PROCESS 


Code Selection Based on: 

• Availability 

• Reliability 

• Efficiency 

Turbulence Model/Near-wall Treatment Selection Based on: 

• Available Database 


Result: 


RNG k-s with Two-laver (Low-Re) Near-wali Model 
and Tetrahedral Cells 
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COMPARISON WITH EXPERIMENTAL DATA (M=0.6 CASE) 


Baseline Configuration (Round 35 Injection) 
MbO.6, U0.23, DR=1.6 



0 5 10 15 20 25 30 

x/d 


Baaollno Configuration (Round 35 Injoction) 
xJM, UoO.6, U0.23, OR»t.6 




Centerline Effectiveness 


Local Effectiveness 






COMPARISON WITH EXPERIMENTAL DATA (M=1U> CASE) 


Baseline Configuration (Round 35 Injection) Baseline Configuration (Round 35 injection) 

Mai .Of 1*0.63, DR=1,6 x/d=10, M»1.0, 1*0.63, DR*1.6 



0 5 10 15 20 25 30 -2 -1.5 -1 -0.5 0 0.5 1 

x/d z/d 


Centerline Effectiveness Local Effectiveness 


COMPARISON BETWEEN ROUND AND NEW SHAPE HOLE RESULTS 


Round vs. NASA/New Shape Results Round vs. NASA/New Shape Results 

11*1.0, 1*0.63, DR-1.6 x/d-10, M*1.0, 1*0.63, DR=1.6 



0 5 10 15 20 25 30 35 

x/d 


Centerline Effectiveness Local Effectiveness 
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Linear Cascade Rig 



TRACER-GAS TECHNIQUE TO MEASURE LOCAL 
FILM COOLING EFFECTIVENESS 


100% co 2 

\ 

* 

I 

To gas analyzer or 
pressure transducer 


Relies on heat-mass transfer analogy 
Heavy gas to match density ratio(up tol .5) 
Sampling taps can also be used as press, taps 
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CLOSURE 


INTERIM RESULTS 

• New hole shape developed which meets program 
goal of increasing TRIT by 1 00 - 200 «F 

• 3-D CFD is an effective design tool for film cooling 
prediction 


FUTURE MILESTONES 

< 

• Develop affordable manufacturing method 

• Experimental validation of new hole shape 
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TURBINE AIRFOIL FILM COOLING— DESIGN INTEGRATION 


Robert Bergholz 

General Electric Aircraft Engines 
Cincinnati, Ohio 


AGENDA 


• Turbine Airfoil Film Cooling - Design Integration 

• GEAE Single-Passage Cascade Facility 

• GE CR&D Transient Cascade Facility 

• Turbine Airfoil Film Cooling Design Issues 


Turbine Airfoil Film Cooling - Design Integration 

Objectives 

• Maximize coolant heat capacity utilization. 

• Focus on overall heat flux management via 
film and thermal barrier technologies to control 
structural thermal gradients. 

• Optimize overall film effectiveness distribution. 

• Customize film hole geometries for improved 
effectiveness with lower blowing ratios. 
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Approach 

• Acquire film effectiveness and external heat 
transfer coefficient data on promising film 
cooling configurations. 

• Apply optimization and computational DOE 
methods to quantify and rank design payoffs. 

• Use unstructured grid CFD tools to predict film 
hole and external flowfields, local heat transfer 
effects, and customized film hole geometries. 


Results 

• GEAE Single-Passage Cascade 

- Gas concentration / mass-transfer method. 

- t| f results for several airfoil geometries. 

• GE CR&D Transient Cascade 

- Both Tj f and h measurement capability. 

- Data acquired for full-scale vane geometries. 

• GE CR&D Water Tunnel 

- PLIF measurement technique. 

- Full-coverage, discrete-hole film cooling 
applicable to combustor liners. 

Milestones 

• 1996 

- Single-Passage Cascade 

• Engine airfoil support tests. 

- Transient Cascade 

• GEAE / CR&D IR&D film tests. 

• 1997 - 99 

- Single-Passage Cascade 

• Engine blade / vane, showerhead, ... 

- Transient Cascade 

• AST film tests. 
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GEAE Single-Passage Cascade 

Fred Buck (GE Aircraft Engines) 


BACKGROUND - FILM COOLING TESTS 

AT GEAE 


• Past History 

- The current airfoil film cooling data base at GE Aircraft Engines 
was obtained using full annular cascades, sector cascades or 
2-D linear cascades run at different mainstream and coolant 
temperatures. 

- For large aircraft engine turbine airfoils (e.g., GE90, CF6-80 and 
CFM56), these tests require large amounts of airflow (e.g., for an 
annular cascade ~ 40 kg/sec) and significant manpower and 
financial resources. 

• Current Reality 

- Resource reductions mandated by the current business climate 
have reduced the available opportunities to run such test 
programs. 


A more economical and efficient test method was 
required to continue to obtain the required r| 
design data base. 
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SINGLE PASSAGE RIG CONCEPT 



Modelling a single passage would significantly 
reduce the mainstream flow required to match 
engine Mach and Reynolds numbers. 


FILM COOLING EFFECTIVENESS 
MEASUREMENTS MADE USING MASS 
TRANSFER/GAS ANALYSIS 


• For thermal film cooling effectiveness 
measurements: 

- Tests are run with mainstream and coolant flows at different 
temperatures. Mainstream and/or coolant temperature should 
be set to match engine coolant-to-mainstream density ratio. 

"fgffcooi " Pcoo/Pg 

- Effectiveness measurements made from wall and flow 
temperatures using 

'H ~ OVec " "^aw.cVOVec ‘ "^cool) 

Wall temperatures in most cases will be measured using 
embedded thermocouples. 

- Corrections may be required to account for nonadiabatic wall 
conditions (e.g., conduction within the cooled wall and/or 
radiation). 
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MASS TRANSFER/GAS ANALYSIS 
TECHNIQUE (CONTINUED) 

• For gas analysis film cooling effectiveness 
measurements: 

- Tests are run using a foreign gas or a mixture of gases as the 
coolant. The engine coolant-to-mainstream density ratio should be 
matched. 

- Film cooling effectiveness measurements are made using foreign 
gas concentrations measured at the wail 

Tl = (^g “ ^wallV(^g ” ^cool) 

- In this test, the wall gas concentrations were measured using static 
pressure taps installed downstream of the film injection site. 
Multiple spanwise tap locations were required to determine the 
average spanwise film cooling effectiveness. 

- The rate at which the wall samples are extracted should be 
controlled. 

Film cooling effectiveness measurements made 
using gas analysis would permit testing with ambient 

mainstream flow . 


SCHEMATIC OF TEST 
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TEST MODEL CROSS-SECTION 



MAINSTREAM PLENUM 

SCALE ( mm > 

0 20 40 60 60 100 120 140 160 180 200 

1 1 1 1 1 1 1 1 I I I 1 I I I I I 1 1 I I 

The test model consisted of 2-D components which 
were sandwiched between two slotted endwalls. 


FILM HOLE GEOMETRY AND TEST 
MODEL INSTRUMENTATION 
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Ps/Pt 


LOCAL PRESSURE DISTRIBUTION - 
COMPARISON WITH DESIGN INTENT 



The cascade pressure distribution was matched 
over the entire single passage datum surface. 


LOCAL RESULTS - PRESSURE SIDE 

m ~ 1.0 



0 10 20304090607080 

x/d 


On pressure side for x/d > 10, In-Line and Offset 
were approximately equal . 
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LOCAL RESULTS - SUCTION SIDE 

m ~ 1.0 



On suction side at x/d ~ 25, In-Line and Offset i\ still 
were measurably different. 


GE CR&D Transient Cascade 

Ron Bunker (GE CR&D) 
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Objective 


• Measure Both Heat Transfer Coefficients and Film 
Cooling Effectiveness Values in Linear Airfoil 
Cascades Under Pressure and Temperature Conditions 
which Simulate Engine Non-Dimensional Parameters 
(Reynolds Numbers, Mach Numbers, Coolant-to-Gas 
Density Ratios, Met Turbulence Intensities) 

• Build and Operate a Transient Test Facility Similar to 
the One Used by Lander, Fish and Suo (J. Aircraft, 
1972) 


Transient Test Methodology 


• Adiabatic Film Effectiveness 



• Heat Transfer Rate (Lander, Fish, Suo, 1972) 

q. = h,(T.»-T.) = plc p ^ 



293 



Transient Test Methodology 


• Curve Fitted Equation from Measured Temperatures and 
Times 

»«rt 

T w (steadystate)-T w (t) =g picp 
T w (steady state) — T w (t = 0) 


Schematic of Test Facility 



PNEUMATIC 

ACTUATOR 
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film EffecHwmess Temperature (C) 


Airfoil Temperature and Film Cooling Effectiveness 

Reynolds Number=5 17,000 (Based on Inlet conditions and Chord 

Length) 


200 
150 
100 
50 
0 

-1 -08 46 <04 -02 0 02 04 06 0.8 1 

Dimensionless Distance from LE. 




-1 -afi * 0.6 - 0.4 - 0.2 0 0.2 04 04 0.8 1 


Dimensionless Distance from LX, 
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Summary and Conclusions 

» Transient Test Facility Built to Determine Both Heat Transfer and 
Film Effectiveness in Linear Airfoil Cascades 

— Steady-state Operation Provides Film Effectiveness 

- Transient Operation Provides Both 

• Base line Tests Conducted with Cylinders In Cross Flow - Stagnation 
Point Heat Transfer Agrees with Established Correlations 

• Linear Airfoil Cascade Tests Show; 

- Mach Number Distributions Agree with Inviscid Code Predictions 
— Airfoil Heat Transfer Coefficients Measured with Plugged Film 

Holes Follow Boundary Layer Code Predictions 

- Film Injection Increases Heat Transfer Levels both on Pressure 
and Suction Sides 


Turbine Airfoil Film Cooling Design issues 

• Expand the current film cooling database to cover a broader range of engine-level 
conditions for both blades and vanes ... 

- Reynolds number 

- Blowing ratio 

- Density ratio 

- Roughness 

• Develop small-scale transducers for direct measurement of heat flux in hlgh- 
temperature environments. 

• Focus on enhanced cooling designs for “edges” - tips, trailing edges. 

• Apply CFD technology to screen film hole geometries. 

- Interior hole flow 

- Film jet vortex - wall interaction 

- Shear-layer / wall-layer dynamics 

• Understand relative film-TBC-convection payoffs for different design requirements. 

- Optimization targets 

• Min Wc at constant life 

• Max SFC (Fn I Weight) at constant part life 

• Max life (min LCF) at constant WC 

• Min life-cycle cost - component / system 

- Inverse cooling configuration design for prescribed mechanical / cooling constraints 

• Create flexible design systems responsive to increasingly rapid engine development 
cycles. 
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AN INITIAL MULTI-DOMAIN MODELING OF AN ACTIVELY COOLED STRUCTURE 


Erlendur Steinthorsson 
ICOMP 

Cleveland, Ohio 

OUTLINE 

• Introduction 

• Objectives of current research 

• Methodology for coolant flow simulations 

• key components 

• key advantages 

• Multiblock/Multigrid scheme - sample results 

• Future research 

• Summary 

INTRODUCTION - The TRAF3D.MB Code 

A computer code designed for simulating flow and heat transfer in turbo- 
machinery: 

• Based on the TRAF3D code (Amone, et al,, 1992), 

• finite volume discretization (second order accurate), 

• multigrid convergence acceleration, 

• multi-block grids for complex geometries, 

• Algebraic and two-equation turbulence models. 
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INTRODUCTION - The TRAF3D.MB Code 


Code has been applied to several turbomachinery flows, including 

• turbine blade and blade tip heat transfer (Ameri, et al., 1995, 1996), 

• study of heat transfer on recessed tips (squealer tips, Ameri, et al., 
1997) 

• study of heat transfer in serpentine passages and channels with ribs 
and bleed holes (Rigby, et al., 1996, 1997), 

• simulations of flow and heat transfer within the “branched duct” 
coolant passage model (Steinthorsson, et al., 1996, 1997) 


Why is the TRAF3D.MB code designed the way it is? 


What is the rationale behind the main code-design choices, 


...to use multi-block grid systems for complicated geometries, 

and 

...to use multigrid convergence acceleration. 
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OBJECTIVES 


• Develop a methodology for 3-D numerical simulations of flow and 
heat transfer in turbomachinery that is 

• capable of delivering accurate results, 

• computationally efficient, 

• capable of handling complicated geometries, 

• efficient in terms of human effort, 

• useful as a platform for evaluating turbulence models, film- 
cooling flow models and other models for turbine cooling 
flows. 

• Demonstrate the effectiveness of the methodology by simulating 
benchmark flows and by computing flow and heat transfer in gas 
turbine engines at realistic operating conditions. 


DEVELOPMENT OF METHODOLOGY 

Key issues: 

• Accuracy, 

• efficiency, 

• modularity. 

Constraint: 

• Must be applicable to flows in complicated geometries. 
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METHODOLOGY FOR 3-D NUMERICAL SIMULATIONS 


Key issues that determine Accuracy and Efficiency: 

Accuracy: 

• discretization of the governing equations 

• quality of grid system 

• resolution of the flow field 
Efficiency: 

• automation of grid generation 

• choice of algorithms 

• optimization of computer codes (vectorization, parallelization) 

• structure of grid system 

• judicious use of grid points. 


TRAF3D.MB 


Key Design Choices: 

• Use multi-block grid systems. 

• Use multi-grid convergence acceleration. 

• Use finite volume discretization. 
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MULTI-BLOCK GRID SYSTEMS 


Globally unstructured: 

• “Unstructured” assembly of blocks, 

• great flexibility for modeling complicated geometries. 

Locally structured: 

• Each block is a “body-fitted” structured grid, 

• grid lines are at least C 1 continuous at block boundaries. 


WHY USE MULTI-BLOCK STRUCTURED MESHES? 

For reasons of accuracy : 

• well suited for viscous boundary layers, 

• high-order schemes can be implemented “easily.” 

For reasons of efficiency: 

• simple array data structures are applicable, 

• high level of optimization on computers, 

• low memory overhead, 

• regular cell-to-cell connectivity allows effective application of 
implicit schemes, multigrid schemes and other convergence accel- 
eration techniques. 
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GENERATION OF MULTI-BLOCK GRIDS 


Current “state of the art”: 

• Commercial grid generator is available. 

• Grid generation is “semi-automatic”; 

• shape of blocks and grid point distribution is optimized, 

• topology needs to be specified. 

• Human effort is greatly reduced compared to older structured grid 
generation technology. 

• Can be fully automated for specific classes of geometries. 


TRAF3D.MB - Ability to Handle Complicated Geometries 

• Completely general multiblock capability: 

• no limitation on number of blocks 

• no limitation on connectivity between blocks 

• Writes blocks to files (SSD) if memory is limited 

• Ability to handle block interfaces with discontinuous grid lines 
across branch cuts. 

• Automatically detects periodic block interfaces 

• linear (e.g., linear cascades) 

• annular (requires rotation of velocity vectors upon transfer across the interface). 
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SAMPLE GEOMETRY HANDLED BY TRAF3D.MB 



MULTIGRID CONVERGENCE ACCELERATION 


• Makes convergence rates essentially independent of grid resolu- 
tion. 

• Reduces required CPU time by an order of magnitude or more. 

• Leads to at most 20% increase in memory 

• Especially useful for simulations requiring high resolution of the 
flow field. 

• Small increase in code complexity (same routines can be used for 
smoothing on all levels). 
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TRAF3D.MB - OVERVIEW 


• Full, compressible Navier-Stokes equations. 

• Baldwin - Lomax and k-co turbulence models. 

• Rotating reference frames (Chima and Yokota, 1990). 

• Discretization (2nd order accurate in space): 

• cell-centered finite volume discretization 

• central differencing (Swanson and Turkel, 1987) 

• central differencing for diffusion terms 

• Runge-Kutta time integration with local time stepping and implicit 
residual smoothing (Jameson, et al. 1981, 1983, 1988), 

• Multigrid scheme (FAS, FMG; Brandt, 1977, Jameson, 1983). 


TRAF3D.MB - Ability to Handle Complicated Geometries 

• Completely general multiblock capability: 

• no limitation on number of blocks 

• no limitation on connectivity between blocks 

• Automatic arrangement of blocks into files if memory is limited. 

• Writes blocks to files (SSD) if memory is limited 

• Ability to handle block interfaces with discontinuous grid lines 
across branch cuts. 

• Automatically detects periodic block interfaces 

• linear (e.g., linear cascades) 

• annular (requires rotation of velocity vectors upon transfer across the interface) 
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TRAF3D.MB Sample Results 
Flow and heat transfer on squealer tips: 

• E3 turbine blade. 

• Results obtained using k-(0 turbulence model. 

Flow and Heat transfer in “branched duct” geometry. 

• Nominal Reynolds number: Re = 335,000. 

• Results obtained using Baldwin-Lomax turbulence model. 

• Results obtained using k-co turbulence model. 


Flow and Heat Transfer on E3 Rotor Blade with Squealer Tip 




Branched Duct - Heat Transfer (Re=335,000) 



(a) 



(b) 


(a) Results obtained with Baldwin-Lomax 
Turbulence model. 

(b) Results obtained with k-co Turbulence 
model. 

(c) Experimental results (NASA TM 
106189, AIAA-93-1797). 



FUTURE RESEARCH 

• Implement additional turbulence models. 

• Add parallel computing capability. 

• Study methods to improve rate of convergence to steady state for 
multi-block grids containing large number of blocks. 

• Implement improved convergence acceleration schemes, e.g, pre- 
conditioning methods. 

• Implement solution adaptive refinement capability. 

• Implement heat conduction module for conjugate heat transfer 
studies. 
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SUMMARY 


• A methodology for the simulation of turbine cooling flows is being 
developed. 

• The methodology seeks to combine numerical techniques that opti- 
mize both accuracy and computational efficiency. 

• Key components of the methodology include the use of multiblock 
grid systems for modeling complex geometries, and multigrid con- 
vergence acceleration for enhancing computational efficiency in 
highly resolved fluid flow simulations. 

• The use of the methodology has been demonstrated in several tur- 
bomachinery flow and heat transfer studies. 

• Ongoing and future work involves implementing additional turbu- 
lence models, improving computational efficiency, adding AMR. 
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AERO-THERMO-STRUCTURAL OPTIMIZATION OF COOLED TURBINE BLADES* 


George Dulikravich 
Pennsylvania State University 
University Park, Pennsylvania 


OBJECTIVES: 

Develop a design system for constrained 
optimization of 3-D cooled turbine blades. 

Provide industry with a modular design 
optimization tool that will take into 
account interaction of the hot gas flow 
field, heat transfer in the blade material 
and the internal coolant flow field, and 
stress/deformation field in the blade . 


EXPECTED BENEFITS: 

1 . fully 3-D analysis and design 
capability instead of 2-D 

or quasi 3-D capability, 

2. simultaneous account of heat transfer, 
aerodynamics, and elasticity 
instead of aerodynamics alone, 

3. maximum use of the existing 
flow-field analysis software, 

4. optimized aerodynamic blade 
shape for minimum total pressure 
loss and maximum torque, 


* Presented at the conference by Joe Gladden, NASA Lewis Research Center. 


311 



5. ability to specify geometric, elasticity, 
flow field and thermal constraints, 

6. optimization of hot surface heat 
flux and temperature distributions, 

7. minimized number of coolant flow 
passages inside the blade, 

8. minimized coolant mass flow rate 
as a result of minimized integrated 
blade surface heat flux, 

9. minimized coolant pressure drop 
due to the resulting smooth wall 
coolant flow passages, 

10. reduced need for manufacture 
of small holes for blade film 
cooling, 

11. reduced need for elaborate 
pattern of pins and fins on 
cooling passage walls, and 

12. minimized design cycle time. 
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PRESENT STATUS: 


We have already developed and 
tested the following software: 

1) Inverse design of single coolant 
passage shapes in 3-D coated 

& uncoated turbine blades, 

2) Inverse determination of 
temperatures, heat fluxes, and 
convective heat transfer coefficients 
on 3-D coolant passage walls or on 
blade hot surface, 


3) 3-D finite element thermo-elasticity 
dynamic analysis of unsteady 
stresses and deformations, 

4) A hybrid constrained optimizer 
incorporating: 

gradient search, 
genetic algorithm, and 
simulated annealing. 

5) Non-Iterative Inverse Determination of 
Temperature-Dependent Thermal 
Conductivities. 

WORK IN PROGRESS: 

1) Optimization of thicknesses of the 
blade wall and multiple internal struts. 
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AERODYNAMICS 

Conic Section Design Variables and Turbine Blade Definition 


Uses conic curves connected between control points 



2D turbine airfoil definition using conic sections and control points 
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Constraint Functions 


inequality and equality constraints 

* 

g k (x)^° h k (x) = 0 

Perimeter 

Cross-sectional Area 
Chord Angle 
Moment of Inertia 
Airfoil Thickness 
Platform 
Leading Edge 
Trailing Edge 
Axial Chord 
Throat Area 
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GENERIC TURBINE ROTOR BLADE 





Aerodynamic shape design of generic turbine rotor airfoil using 
conic section design variables and target pressure distribution function. 
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SHAPE OPTIMIZATION OF 
3-D COOLANT FLOW PASSAGES 



3 


L2N0RM (%) 



Convergence history of the objective function for the 
three-dimensional turbine blade wall thickness optimization. 




HEAT TRANSFER 


WHAT IS INVERSE HEAT CONDUCTION? 
DETERMINE 

Temperatures, Heat Fluxes and Convective Heat Transfer 
Coefficients on Unknown, Inaccessible Surfaces 

GIVEN OVER-SPECIFIED INFORMATION 

Interior Temperature and Temperature Gradient Measurements 

Over-Specified Thermal Boundary Conditions on Accessible Surfaces 

(Temperature and Heat Flux) 

or Q = hlT-T^b) 
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APPROACHES 


ITERATIVE 

Sum of Least Squares Minimization 
Errors are Magnified 

Requires Regularization and Smoothing of Extrapolated Heat Fluxes 

DIRECT (IMPLICIT) 

Forms a Singular Solution Matrix 
Discard Eigenvalues That Amplify Errors 
Uses SVD 


320 



APPLY BOUNDARY CONDITIONS 


Example 



u = Uj <& q =q 2 u = u 2 <& q = q 2 


BEM Solution Set 
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Collected Unknowns on the Right-Hand Side 
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Near-Singular or Ill-Conditioned Systems of Linear Algebraic Equations 

[A]{X} = {F} 
















Annular Disk With Interior Temperature Measurements 

? 




Non-dimensional 
temperature 
7 0054344 

6 0043368 

5 00223428 

4 IE-12 

3 -04)123143 

2 -04)296428 

1 -04)469714 


Non-dimensional 

temperature 

7 00593573 

6 00485429 

5 00367857 

4 04)250286 

3 04)132714 

2 IE-12 

1 -04)102429 



Non-dimensional 

temperature 

7 04)525 

6 04)431 

5 04)337 

4 04)243 

3 04)149 

2 04)055 

1 IE-12 


76 Equations 
108 Unknowns 


78 Equations 
108 Unknowns 


81 Equations 
108 Unknowns 


Isotherms Predicted With the Inverse BEM When Temperatures Were Specified at 

4, 6 and 9 Equidistantly Spaced Interior Points 
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Cube-Withln-A-Cube 




Q 

Isotherms on a Slice Through Center of the Cube Predicted by the Inverse BEM 
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Three-Dimensional Rocket Nozzle Wall Section 



Surface Grid of a Rocket Nozzle Wall Section Composed of Four Materials and 

Having a Finned Coolant Passage 



Zirconium 

Oxide 

Nickel- 

chromium 

Copper 

Copper 

closeout 

8.0 

23.0 

378.0 

385.0 


Table. Thermal Conductivities of Different Material Layers 


Tshroud — 140 K 
Ota = -lOOxlO 6 W/nr 
Tamb = 50 K 
hcoav = 1x10 s W/m 2 K 
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Level 0 

F 1923.5 

E 826.557 

0 746.781 

C 669.452 
B 552.454 
A 491.619 
9 437.355 

8 384.202 

7 326.679 

6 306.872 

5 276.419 

4 251.383 

3 217.076 

2 177.857 

1 141.466 


Level 0 

F 1923.5 

E 826.567 

D 746.781 

C 669.452 

B 552.454 

A 491.619 

9 437.355 

8 384.202 

7 326.679 

6 306.872 

5 276.419 

4 251.383 

3 217.076 

2 177.857 

1 141.466 


Isotherms Predicted by the Well-Posed and Ill-Posed BEM for a Three-Dimensional 
Segment of a Rocket Thrust Chamber Wall with a Finned Coolant Passage 
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hconv = 3000 W/m 2 K 
Tgmb — 300 K 
g = 100,000 W/m 3 
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Isotherms (Solid Lines) Surface Grid (Dashed Lines) 


Level t 

8 301,5 

7 301.35 

6 301 2. 

5 301.05 

4 300.9 

3 300.75 

2 300.6 

1 300.45 



Isotherms on the Finned Heat Sink Isotherms on the Finned Heat Sink 
Predicted by Forward BEM Predicted by the Inverse BEM 




Level t 

8 382.073 

7 371.996 

6 36132 

5 351.843 

4 341.767 

3 331.69 

2 321.614 

1 311-537 


Level I 



381.451 

371.452 

361.453 

351.454 
341.456 



Isotherms on the Heat Spreader Isotherms on the Heat Spreader 

Predicted by the Forward BEM Predicted by the Inverse BEM 
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Isotherms (Solid Lines) Surface Grid (Dashed Lines) 


Level t 
8 
7 
6 
5 
4 
3 
2 
1 


381.278 

376.067 

370.856 

365.644 

360.433 

355222 

350.011 


Level t 
8 387.068 

381.696 
378.323 
370.951 
365.579 
360207 
354234 
349.462 


7 

6 

5 

4 

3 

2 

1 






Isotherms on the Substrate Plate 
Predicted by Forward 


Level t 
8 387.068 


Isotherms oil the Substrate Plate 
Predicted by the Inverse BEM 


Level t 
8 386.489 




Isotherms on the Bottom Surface 
Predicted by the Forward BEM 


Isotherms on the Bottom Surface 
Predicted by the Inverse BEM 
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Level t 

8 392.127 

7 391.814 

6 391.5 

5 391.187 

4 390.873 

3 390.56 

2 390.246 

1 389.933 



Level t 

8 391.399 

7 391.098 

6 390.797 

5 390.496 

4 390.194 

3 389.893 

2 389.592 

1 389.291 



Isotherms on the Surface of Chip Isotherms on the Surface of Chip 

Predicted by the Forward BEM Predicted by the Inverse BEM 
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Conjugate Thermal Entry Problem 

y=+h v x =0 T=Ti 



Region 1 : Solid with conductivity k s 
Region 2: Fluid with conductivity k F and known velocity profile v x (y) 
Region 3 : Solid with conductivity k s 



F 1 

E 0.9285715 

D 0.8571429 

C 0.7857143 

B 0.7142857 

A 0.6428572 

9 0.5714286 

8 0.5 

7 0.4285715 

6 0.3571429 

5 0.2857143 

4 0.2142857 

3 0.1428571 
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Figure 3: Transverse dispiacexnents(w) in meters. 


Figure 4: Axial stress(xx) in Pa. 


Level xx 
F 1.0S56E8 
E 8.94Q6E7 
D 9224BB7 
C 8-5091 S7 

S 7.7S33E7 
A 7.0776E7 
9 6.3618E7 

8 5*6461 E7 

7 4.S304E7 

S 4J2146E7 
5 3.4989 E7 

4 2.7831 E7 

3 2.0874E7 

2 1.3517E7 

1 S*3S9S€6 


-6.937E-6 
C -3J25E-6 
B -1JS6E-5 
A *1 .387E-5 
9 -1.61SE-5 

8 AJB9&6 
7 -2.081 &6 

6 -2-312E-5 

5 -2S43E-5 

4 -2.77SE-5 

3 -3.Q06E-5 

2 -3.237E-5 

1 -3.468S-5 







GROUP DISCUSSION - December 12, 1996 

SESSION IIL COOLANT FLOW INTEGRATION AND OPTIMIZATION 


Rob Simoneau - Facilitator Jim Heidmann - Scribe 

Nirm Nirmalan to give review of work at Allison. 

Peer Review by industry told NASA to take a leadership in interdisciplinary/ 
multidisciplinary design analysis methods - a system that couples the internal cooling 
system through the film cooling holes to the external aerodynamics (1993). 

Industry was telling NASA to “stay out of’ industry design systems, but to provide tools. 

Early work used quasi, 3D, invicid programs (TSONIC) on the external side and network 
nodal analysis with correlations on the inside. 

Bob Simoneau knows of no system to tie the whole problem together. 

Joel Wagner 

People still use meanlines and streamlines to correlate results. There is a blurring of the lines 
between primary and secondary flow paths. He doesn’t see 3D driving the design, but inputting 
information. 

Rob Norton 


Need to look at effect of variations in shapes of pedestals, pins, etc. to study sensitivity to small 
changes in manufacturing tolerances. Look at trade-offs. 

Nirm Nirmalan 


Alli son - cost modeling of airfoils. Designer should be able to estimate cost of design early in 
design process. Most cost is decided in early design - detailed design has only a small influence. 
Would like to put this info into design codes. Get data from experimental cascades, etc. Getting 
heat transfer designer to think about manufacturing (especially casting) is difficult. 

Bob Simoneau 


Asked if thermal design is equally important to aero in design? 
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Bob Bereholz 

Seeing less emphasis on narrow specialization, particularly in gas path - aero designer considers 
heat load. “Aerothermal” is becoming unified. More of this in future. 

Rob Norton 

Agrees with above point. Need to look strongly at manufacturing because that is where the big 
costs are. 

Boris Glezer 

“Inter-discipline” is key word. Small companies learned this early. EPS showing more signs of 
interaction in industry: e.g. endwall cooling, combustion interaction, materials, and cooling 
penalties (discharging air requires foil engine testing). Working with universities. 

Bob Bergholz 

Manufacturing - more tools that can demonstrate real payoff - “Can we build it?” Quantitative 
answers. Casting technology has improved through identified payoffs. More willing to invest in 
a process if there is an identified payoff. 

Nirm Nirmalan 

Manufacturing is very uncertain process. Payoff - difficult to identify low NO x payoff. 

Ed North 


Dollar value on NO x - it’s a “go - no go” value. No compromises! It can’t meet NO x standards, 
you can’t fly. Limited by flame temperature. 

Bob Simoneau 


Steady increase in industry use of CFD for design decisions. We are pushing the envelope, lower 
tolerances - is that the right direction? What about fidelity and resolution on the big scale? 

Ron Bunker 

When marketing a power turbine with a combined cycle efficiency of 60%, reaching the goal 
efficiency is highly important. 

Bob Simoneau 


So high fidelity resolution on the big scale is important, right? 
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Ron Bunker 


Yes, integrating aerodynamics, heat transfer, performance, mechanical, life, materials, and 
manufacturing are all important; especially, in steam power turbines. 

Rob Norton 


Question for Yong Kim - Did you compare heat transfer predictions? Answer - No, because of 
tetrahedral elements, but I’m waiting for future improvements in resolution. 

Bob Simoneau 


Have we run out of ideas? 

Bud Lakshminaravana 

New data coming from Europe - detailed measurements of FC effectiveness. We should utilize 
this data. Need to simulate Re, M, Ro, curvature, etc. All have effect on film cooling. CFD is 
bringing Aero and HT committees together. Aerodynamics, manufacturing, and life coming 
together have a long way to go! 3D, Coriolis effects need to be looked at. Need coupled 
plenum, hole, external flow computations to develop models. More experiments in rotating rig. 
What about Wisler’s rig at GE? It’s not being used for film cooling. 

Dave Rigbv 

Why don’t we look at slots or wide holes approaching slots? Less lift-off. 

Boris Glezer 


Structural problem - trailing edge to approach this geometry - it may be possible on pressure or 
suction side. Easy if a simple shape. 

Luzeng Zhang 

Limited cooling air - constrained in spanwise direction. Round hole acts as a meter - a slot 
would not. Laser cannot drill slot. 

Nirm Nirmalan 


Manufacturing process matters - Allison casts holes - slots may be better in a cast blade, but it 
depends on the initial manufacturing process, which is decided by the Aero preliminary design. 
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Ed North 


Endwall cooling - don’t have high pattern factors in steam turbines - endwails are difficult to 
cool. Need cooperation with aircraft engine designer. Can’t afford empirical endwall cooling - 
need to be smart. 

Boris Glezer 


On board: 

Heat Transfer occupies a region inside 
a triangle whose vertices are: 

$/HP, Efficiency, and Life. 

“Interdisciplinary” 

J. C. Han 


Endwall cooling - would linear cascade be sufficient? (to Ed North) 
Ed North 


Need to talk to aerodynamicist. 

Luzeng Zhang 

Can get some info, with 2D cascade since endwall is 3D, but not all. 

Bud Lakshminaravana 

RAE has some data on endwall cooling - cascade must have right velocity profile entering - can’t 
simulate rotor endwall flow in cascade due to rotation, but can simulate stator - Bob B. agrees. 
Does industry have concepts which will do away with film cooling? 

Ed North 

Closed loop cooling is coming for industrial turbines (for ground power) along with combined 
cycle efficiency improvements. Thermal barrier coatings need to be viable for long periods of 
time. CLC is secret to high performance - no film cooling - difficult, but rewards are high. 

Luzeng Zhang 

We aren’t “looking for” film cooling - sometimes we “need” film cooling. They are necessities. 
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Bob Bergholz 

Heat flux reduction is critical. Don’t look at max. temp. - gradients are important. FC controls 
heat flux in blade. A high performance FC design may be bad due to thermal gradients. How to 
control heat flux in part. 

Bob Simoneau 


Good points on what is important in experiments - not much discussion on integration. Why? 

J. C. Han 

It is important, but very difficult to specify appropriate external conditions. Need “real” turbine 
passage for rotating film cooling. These are very difficult experiments, but you need integration 
between the FC and inside supply flow for reliable results. 

Bob Simoneau 

We agree that it is important, then. 

Yung Kim 

Problem is formidable. Codes require a lot of issues - network problem - don’t have good 
database for film cooling discharge coef., need each problem to be subroutine for integration of 
parts. 

Joe Gladden 


ASU and PSU are working on optimization programs. We may optimize to zero coolant. 

Rav Gaugler 

Thanks to all of you for participating in this Workshop and to OAI (especially Karen Balog) for 
hosting it. 

Adjourn. 
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